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PURPOSE 


One  exanple  of  a  SONAR  equation,  for  the  passive  detection  of  a  target 

which  is  radiating  sound  into  the  ocean,  is  ' 

II 

1 

1 

TL  ^  NSL  +  10  log  w  -  DI  +  DT  i 

Source  Level  of  the  target  being  detected  passively 
Transmission  Loss  as  the  signal  propagates  to  the  detector 
Noise  Spectrum  Level  of  the  ambient  noise  in  the  ocean 
Frequency  Bandwidth  of  the  detecting  system 
Directivity  Index  of  the  detecting  system 
Detection  Threshold  of  the  detecting  system  for  a  specified 
probability  of  detection  and  a  specified  false-alarm  rate 

The  purpose  of  this  course  is  to  provide  the  physical  and  conceptual 
background  necessary  to  understand  the  meaning  of  each  of  these  terms 
and  how  they  are  combined  to  form  the  SONAR  equations.  At  the 
end  of  this  course,  the  student  should  be  able  to  apply  the  appropriate 
SONAR  equation  to  a  given  problem  and  estimate  such  things  as  mximum 
detection  range,  optimum  receiver  depth,  necessary  bandwidth,  etc. 


EL  - 
where 
EL  = 
TL  = 
NSL  = 
w  = 
DI  = 
DT  = 
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JotTi  I  Pk  e  s  s  u  r  e  <p(i) 


ACOUSTIC  PFcEGSURE 


Acpc^iiic  t  ’  (P(0  ”  --f(i)  ■ 

f(-i)  -(he  t|^wa  i-l:  y  scytseJ 

ivi  e>$'t  h  'y  d  ^  0  pI'  e>y)€  ^ 


HYDROPHONE  OUTPUT  VOLTAGE 


Hydrophone  sensitivity  'x  =  ^  Ci) 

where  ig  the  instantaneous  voltage  at  output 

terminal  of  hydrophone 


h 


hydrophone 


DISPLAYS 


OSCILLOSCOPE 

instantaneous  voltage  v(t) 


VOLTMETER 

some  average  voltage 

or 
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BEARING-RANGE  RECORDER  (active  sonar) 

brightness  related  to  the  strength  of 
the  echo  returned  from  each  range  and 
bearing 


LOFAR-GRT^  (passive  sonar) 

paper  is  electrically  burnt,  with 
darkness  related  to  the  strength  of 
the  signal  received  at  each  frequency 
cind  time 


MONOFREQUENCY  SIGNALS 


For  ease  of  discussion,  we  will  at  first  restrict  ourselves  to  consideration 
of  tonal s  (monofrequency  signals) .  Later  on,  we  will  see  that  more  compli¬ 
cated  signals  can  be  broken  down  by  FFT,  Fourier  Analysis,  or  other  filtering 
techniques  into  a  collection  of  monofrequency  signals. 

The  time  history  of  the  acoustic  pressure  for  a  monofrequency  signal  can  be 
represented  as 

jc  -■ 

A 

F 

-  1^ 


P  =  The  effective  amplitude  of  the  acoustic  pressure.  This  is  not  the 
same  as  the  peak  amplitude  A,  but  is  more  convenient.  We  will  never 
use  the  peak  amplitude,  but  it  turns  out  that  P  =  0.707A  for  the  above 
signal. 

T  =  The  period  of  the  mono frequency  signal.  The  period  is  the  time 

interval  required  for  the  acoustical  signal  to  go  through  one  complete 
cycle,  returning  to  the  configuration  it  had  at  the  beginning  of  the 
interval . 


T 


T/Al£ 
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PLANE  WAVES 


jc  ~  -S/y  i  -^  ] 

If  X  =  constant  (For  simplicity,  take  x  =  0.) 

A 

/\  - - 


If  t  =  constant  (For  simplicity,  take  t  =  0.) 


I 


As  time  increases,  the  wave  moves  towards  larger  x. 
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it  9  ^ 


A  WAVS  II  SPACS  AND  TIM 
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INTENSITY 


'  -7 


Intensity  =  power  per  unit  area 

where,  "area"  is  the  area  of  a  "window"  perpendicular 
to  the  direction  of  travel  of  the  wave. 

In  terms  of  measureable  quantities, 


where  P  C 

J 


pc 

PC 


=  specific  acoustic  impedance  of  the  fluid 


1.5  X  10  Rayl 
415  Rayl 


for  water 
for  air 


1  Rayl  =  1  kg/m  sec 
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SIZES  OF  THINGS  IN  WATER 


Speed  of  sound,  frequency,  and  wavelength: 


The  speed  of  sound  in  sea  water  is  approximated  by 

3 

c  ~  1. 5  X  10  m/sec  a  few  per  cent. 

It  must  be  remembered  that  because  of  depth,  temperature,  salinity,  and 
contamination  there  may  be  variations  up  to  a  few  per  cent  in  the  exact 
value  of  c. 


Approximate  sizes  of  wavelengtirs  of  monofrequency  sounds  of  different 


frequencies  are: 


Frequency 
10  Hz  =  10 


cps 


50 

100 

500 


1  kHz  =  1000  cps 
5 

10 

50 

1  m  =  1.09  yd 


Wavelength 
150  m 
30 
15 
3 

1.5 

0.3 

0.15 

0.03  =  3  cm 


Pressure  amplitudes: 


2 

1  atmosphere  =14.7  lb/ in  at  the  ocean  surface 
0 

10  yb  (yb  =  microbar) 


(? 


o 


=  lO^^yPa  (yPa  =  microPascal) 
increases  by  about  1  atm  for  every  33  ft  of  depth  in  water. 


Representative  acoustic  pressures  generated  near  a  sonar  transmitter: 

P  ^  10^  to  10^  yb  or  0.01  to  1  atm. 

If  P  --  then  cavitation  is  likely  and  perfomance  nay  be  degraded. 

Representative  weak  acoustic  pressures  detectable  by  a  sonar  receiver: 

P  10”^  yb  or  10~^  atm. 
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SOUND  PRESSURE  LEVEL 


If  P 


If  P 


If  P 


SPL  =  20  log(P/P^^)  re  P^^ 


0.0002  yb  (old  noise  measurements) 

P^^  =  lyb  (old  conventional  standard) 

lyPa  (new  Navy  standard) 


=  1  atm  =  10 
r 

194 


=  10~*^  atm  = 

r 

114 

SPL  =<40 
^140 

=  10"^  atm  = 


t34 
-40 
60 


Exanples 

^yb,  then 
dB  re  0.0002yb 
dB  re  lyb 
dB  re  lyPa 

10  ^yb,  then 
dB  re  0.0002yb 
dB  re  lyb 
dB  re  lyPa 

10”^yb,  then 
dB  re  0.0002yb 
dB  re  lyb 
dB  re  lyPa 
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CONVERSION  BETWEEN  SPL’s  WITH  DIFFERENT  REFERENCE  PRESSURES 


For  a  given  pressure  amplitude 


SPL  re  lyPa  =  SPL  re  lyb  +  100  dB 
SPL  re  0.0002yb  =  SPL  re  lyb  +  74  dB 


Notice 


0  dB  does  not  mean  the  absence  of  sound  ,  but  neans  P  =  P 

-  ref 


If 

then 


P  =  P 


ref  ' 


SPL  =  20  log  (P/P^f)  =  20  log(P^/P^f)  =  20  log  1  =  0  dB  re  P^^ 
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TKANSmSSION  LOSS 


PCr^) 


SPL(rj^) 

=  20  log 

ref 

SPL(r2) 

P(r.) 

=  20  log  p-^ 

ref 

SPL(r^) 

-  SPL(r2)  “ 

PCr^^) 

P(r2) 

TL  =  SPL 

(1  meter)  -  SPL  (r  meters)  = 

P(l) 

KrT 
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SAMPLE  CALCULATION  OF  TRANSMISSION  LOSS 


If  SPL  at  1  meter  = 

120  dB  re  lyb 

and  SPL  at  range  r 

=  40  dB  re  lyb 

then  TL  =120  -  40  dB  =  80  dB 

TL  does  not  depend  on  the  reference  pressure.  Let's  convert  the  above  SPL's 
to  reference  lyPa  and  recalculate  the  TL 


New  SPL  at  1  meter  = 

220  dB  re  lyPa 

and  SPL  at  range  r  = 

140  dB  re  lyPa 

so  that  TL  =  220  - 

140  =  80  dB 

14 


TRANSMISSION  LOSS  FROM  SPREADING 


The  total  transmission  loss  TL  can  be  considered  to  arise  from  two  different 
causes : 

1)  As  the  rays  of  sound  propagate  out  from  the  source  and  travel  throu^ 
the  water,  they  will  bunch  together  or  spread  out  depending  on  the  properties 
of  the  speed  of  sound  profile  for  the  water.  The  transmission  loss  which  arises 
from  this  effect  is  termed 

TL  =  transmission  loss  from  spreading 
Some  sinple  examples  of  spreading  for  rays  traveling  straight  lines  are 

Spherical  spreading 

TL  =  20  log  r 
g  ^ 


Cylindrical  spreading 
TL  =  10  log  r 

s 


A 

No  spreading 

TL  =  0  dB 
g 


In  the  real  ocean,  the  rays  never  travel  in  strai^t  lines,  so  that  these 
above  sinple  examples  are  unrealistic.  We  will  see,  however,  that  there  are 
many  situations  in  which  the  sound  will  spread  out  spherically  or  cylindrically 
over  large  portions  of  its  path,  and  the  above  simple  equations  will  appear 
in  the  complete  transmission  loss  equation. 
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SPHERICAL  AND  CYLINDRICAL  WAVES 


SPHERICAL 


CYLINDRICAL 


jb  -  ”  SI  A/  (i  -  ^ 

IF  'T  »  >v 
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TRANSMISSION  LOSS  FROM  DISSIPATION 
2)  There  are  a  variety  of  other  physical  effects  termed  dissipative 
effects,  vhich  reduce  the  SPL  of  the  sounds  in  addition  to  the  spreading 
of  tl^  rays.  Since  the  spreading  of  rays  is  described  by  TL  ,  we  will  call 
the  portion  of  the  transmission  loss  which  results  from  these  other  effects 
TL  =  Transmission  Loss  from  dissipation 

X/ 

There  are  a  number  of  these  which  are  important  in  underwater  sound  propagation: 

(a)  Absorption.  As  sound  tra.vels  throu^  the  water,  there  is  a 
conversion  of  the  acoustical  energy  into  thermal  energy  because 
of  effects  vhich  are  very  similar  to  friction  in  moving  machinery 
and  hysteresis  losses  in  transformers. 

(b)  Transmission  into  the  bottom:  When  a  sound  beam  strikes  the  ocean 
bottom,  part  of  the  beam  is  reflected  back  upward  in  the  water 
and  another  part  is  sent  dcwn  into  the  bottom.  The  beam  vhich 

is  reflected  back  upward  in  the  water  contains  less  energy 
than  the  beam  which  traveled  down  to  the  bottom. 

(c)  Scattering.  The  presence  of  rough  boundaries  at  the  ocean  surface 
and  bottom  can  lead  to  a  scattering  of  the  sound  out  of  the  reflected 
beam.  The  presence  of  bubbles  and  other  impurities  in  the  water 
itself  can  also  scatter  sound  energy  out  of  the  path  of  the  ray. 

(d)  Leakage.  In  certain  kinds  of  sound  propagation  in  the  ocean, 
like  the  surface  layer,  the  rays  of  sound  can  "leak"  out  the 
bottom  of  the  layer  each  time  they  bend  down  to  graze  it.  This 
"leakage"  results  in  a  loss  of  energy  from  the  layer. 
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TRANSMISSION  LOSS 


In  general  the  total  transmission  loss  TL  can  be  composed  of  tvro  terms. 


TL  =  TL 

+  TL. 

g 

1 

where 


TL  =  transmission  loss  from  geometrical  spreading 


TLjj^  =  transmission  loss  from  dissipation 
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These  curves  show  how  the  speed  of  sound  is  related  to  depth  and 
temperature.  The  effect  of  salinity  is  not  included,  so  these 
curves  are  appropriate  only  for  vjater  of  constant  salirity, 
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SPEED  OF  SOUND  POEMUIA 


C  I 


-t- 

-f 

1.3  (5-3?) 

-f 

o,o\n  y 


Temperature  Correction 
T  in 


Salinity  Correction 
S  in  parts  per  thousand 


Depth  Correction 
z  in  meters 


Very  accurate  formulas  for  the  speed  of  sound  as  a  function 
of  temperature,  salinity,  and  depth  are  very  involved  and 
difficult  to  calculate.  The  above  expression,  although 
not  highly  accurate,  is  sufficient  for  many  calculations 
and  is  relatively  easy  to  use.  It  is  accurate  to  within 
about  0.1  meter  per  second. 

For  a  discussion  of  many  of  the  more  precise  formulas  cuid 
the  discrepancies  between  them,  the  interested  reader  is 
referred  to  V.  A.  Del  Grosso,  J.  Acoust.  Soc.  Am.56 , 
1084-1091  (1974) . 
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TYPXA'^  soum  SPEW  PROFILE 
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PROFILE  VAPJATIOIS 
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Depth 


Speed  of  Sound  c(z) 


RAIS 


Ran^e  r 


RATE  IN  CONSTAOT  GRADIENTS 


c(^) 

C. 


f 


is  the  £,radient 


The  gradient  g  can  be  either  positive  or  negative. 

For  a  speed  of  sound  profile  composed  of  layers  vathin  each  of 
vhich  the  gradient  has  constant  value,  then  the  segment  of  the 
ray  'V'ithin  each  layer  is  the  arc  of  a  circle  whose  center  lies 
at  the  depth  z  ■  -c^/g  below  the  top  of  the  layer  (x^here  the 
speed  of  soxxnd  is  c^  at  the  top  of  the  layer  and  g  is  the  value 
of  the  gradient  for  that  layer).  The  radius  of  the  circle  is 

R  '*  I  ^  - 

^  ~  I  & 

The  ray  always  bends  toward  the  lower  speed  of  soxind. 
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CONSTRUCTION  OF  RAY  PATffi 


Line  of  centers  for  bottom  layer 


I 


I 

! 

I 

! 
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ACOUSTICAL  RECIPROCITY 


If  the  direction  the  sound  travels  over  the  propagation  paths  is 
reversed  by  exchanging  the  positions  of  the  source  and  receiver,  and  the 
propagation  conditions  remain  the  same,  then  the  transmission  loss 
between  source  and  receiver  will  remain  the  same. 


Active 


As  a  consequence,  the  transmission  loss  suffered  by  the  signal  in 
reac'ning  the  target  from  an  active  source  will  be  identical  with  the 
transmission  loss  experienced  by  a  signal  of  the  same  frequency  (and 
duration)  radiated  from  the  target  when  it  reaches  a  receiver. 
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RAYS  AND  A  SUTtFACE  DUCT 


LOSS  FROM  THE  SURFACE  DUCT 


When  a  ray  reflects  from  the  ocean  surface,  some  of  the  energy 
is  scattered  into  different  directions  and  can  propagate  out  of  the 
surface  layer.  This  energy  is  thus  lost  from  the  channel,  which 
increases  the  transmission  loss.  Additionally,  energy  can  diffuse 
out  of  the  bottom  of  the  layer.  These  effects  can  be  described  in 
terms  of  the  "loss  per  bounce".  The  number  of  bounces  is  the  range 
divided  by  the  skip  distance  r^. 

An  empirical  formula,  resulting  from  the  analysis  of  a  number 
of  experiments,  has  been  obtained  by  W.  F.  Baker,  J.  Acoust.  Soc. 
am.  57,  1198--1200  (1975).“ 


b 


b  =  loss  per  bounce  in  dB 


SS  =  Sea  State 


This  equation  is  valid  within  a  restricted  range  of  parameters, 
given  approximately  by 


2  <  SS  <  5 


3  kHz  <  f  <  8  kHz 


80  ft  K  layer  depth  K  200  ft 


30  ft  <  source  depth  <  60  ft 
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TRAffiMISniON  LOSS  IN  THE  KLXED  UYER 


The  depth  of  the  mixed  layer  is  designated  by  D. 


The  depth  of  the  source  in  the  layer  is  designated  by  z  . 

5 

The  transition  range  r  is  the  distance  the  sovind  beam  must  travel 
until  it  effectively  ^  fills  the  channel.  It  is  calculated  from 


A  -  {.05'*I0 
t 


^  I  ^ 
V  ■D-lb 


(all  distances  are  in  meters) 


The  skip  distance  r^  is  the  distance  required  for  the  rays  vhich  just 
graze  the  bottom  of  the  layer  to  complete  one  full  cycle  from 

the  layer  bottom  to  the  surface  and  back  to  the  bottom.  It  is 
calculated  from 

A-  -  (distances  are  in  meters) 

A  ^ 


If  the  loss  per  bounce  is  designated  by  b  and  the  energy  loss  resulting 
from  absorption  is  designated  by  a,  then  the  transmission  loss  for  the 
mixed  layer  can  be  estimated  by  tTie  folloi'.dng  formulas: 


TL  ■  20  log  r  +  ar  within  transition  range 

TL  ■  10  log  r^  +  10  log  r  +  ar  +  tir/r^  beyond  transition  range 


For  ranges  less  than  r  the  sound  is  spreading  out  spherically  to  fill 
the  sound  channel.  For  ranges  greater  than  r  ,  however,  the  sound 

has  spread  far  enough  vertically  to  be  reflected  ^  from  the  surface  and 
bent  back  upward  from  the  bottom  of  the  layer.  Once  this  hap  ;ens,  then 
the  sound  is  trapped  and  will  spread  cylindrically. 
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C4LCUL/i'jI0N  OF  A 


TRAFFinSSIOI,  LOSS 


Assupiptioiis: 


Depth  of  the  traxed  layer  =  itO  in 

Deptli  of  the  source  *  20  n 

Frecuency  of  tl'O  sound  ■»  2  kHz 

Sea  state  «=  3 

1  e  then  determine 

Skip  distance  =  ^300  m 

Transition  ranpe  =  9u0  m 

Attenuation  constant  =  0.00015  dE/rn 

Loss  per  bounce  -  3.5  dB/bounce 

Tlie  transmission  loss  formulas  are  thus 

(a)  For  ranj^es  less  tlian  the  transition  ran^e, 

TL  -  20  loL  r  +  (l,5  x  10"^)r 

(b)  For  rar£es  {renter  than  the  transition  ran^e, 

TL  =  10  lop  r  +  (8.1  X  10‘’^)r  +  30  . 
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TRANEMIESIOli  LOSS  CURVE  FOR  SOURCE  AIJD  RECEIVER  IK  THE  DUCT 


Parameters  j 

Depth  of  the  mixed  layer  ^  UO  m 

Depth  of  the  source  -  20  m 

Frequency  of  the  sound  »  ?  kHz 

Sea  state  *  3 

Receiver  is  in  the  duct 


tao  looo  to.oao 


.  (nr 


BUT:  These  results  can  be  radically  chanj^ed  by  effects  such  as 
Surface  Interference  (vhich  becomes  important  at  lower 
frequencies  and  lower  sea  states),  steep  gradients  below 
the  layer,  and  source  and/or  receiver  near  the  top  or 
bottom  of  the  duct.  Additionally,  if  the  frequency  of 
the  sound  is  too  low,  it  will  not  be  trapped  in  the  duct. 

As  a  very  rough  estimate,  the  frequency  must  be  greater  than  about 
(2  X  105)/d3/2,  where  D  is  in  meters  and  the  frequency  is  in  Hz, 
for  the  sound  to  be  well  trapped  in  the  layer  and  for  the  above 
model  to  apply. 
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surface  IhlTERFERENCE 


m 


soi/acE 


RtCttVif^ 


I 


r 


> 


h  ~  SOL'R'C:  Oi^pjft 

d  -  jRECUiiysP  Dtf  '-V 

f  -  {  '-/oRi^orjTAi) 
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Acoustic  pressure  at  the  receiver 


pressure  amplitude  *  A 


>  Ar 


Ad  :3:A 

Ad  »  A 


r 

r 


because  (1)  reflected  rays  travels  a  lon£er  path 
(2)  scattering  from  rough  surface 
if  range  large  and  surface  smooth 
if  surface  very  rough 


phase  delay  ■  ^  (due  to  difference  in  path  length  and  to  phase  change 

upon  reflection) 
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Iransmission  loss 


(-  f  ^  smooth  surface  7^  A^. 


ran£;e  (km) 


range  (Ion) 


rteal  sound  speed  profiles  mil  affect  TL  only  at  large 
distances,  (lotice  the  shadoi':  zone  predicted  at  d  km.) 
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The  rar.g,e  at  vjhich  the  received  si^-nal  is  a  minimum  is  {.iven  by 


r 


2.  hd  c 

me  * 


vhere  'll  2  y  2 ,  '  '  ' 


It  is  possible  to  use  inforuiation  from  the  received  signal  to  determine  the 
depth,  speed,  and  distance  of  closest  approach  of  the  source,  f’or  a  source 
moving  at  constant  speed  v  along  a  straight  line  at  a  fixed  depth  h,  the 


v’here  at  t  =  0,  r  ■  s  ■  the  distance  of  closest  approach. 

Combining  these  equations  gives  the  time  at  which  tlie  signal  is  a  mininum  as 
function  of  frequency. 


a 
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The  display  that  results  is  shown  below: 
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Our  equation  can  be  rewritten  in  standard  form 


where 


O  r 


^  c  S 


I 


and  ^  “  '^/'Zr 


This  is  the  equation  of  a  hyperbola 


where  a  is  the  value  of  f  at  t 


0  and  t/a  is  the  slope  of  the  asymptote. 


If  V  is  calculated  from  the  measured  Doppler  sydfts  of  the  discrete  lines,  * 
the  depth  of  the  source  is  found  from 


and  the  distance  of  closest  approach  from 


Of  course,  for  day-to-day  use,  nomograms  can  be  developed  to  simplify  the 
calculations . 


*  It  will  be  shown  later  that 

where  JLj  is  the  total  frequency  change  of  a  line. 
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Depth 


PROPAGATIOIJ  PATHS  FOR  SHALLOW''  SOURCE  Ai\D  RECEIVER 


(a)  Surfade  Duct.  The  isothermal  water  at  the  ocean  surface  forms 
the  mixed  layer  in  vjhich  sound  can  be  trapped.  The  continual 
reflection  at  the  rou£h  surface,  refraction  at  the  bottom  of 
the  layer,  and  leaka£e  provide  relatively  large  energy  losses 
so  that  this  path  is  useful  for  relatively  short  ranges. 

(b)  Convergence  Zone.  Vhen  the  water  is  deep  enough,  the  sound  rays 
which  travel  down  past  the  SOFAR  channel  axis  can  be  bent  back 
upward  before  they  reach  the  bottom.  These  rays,  when  they 
approach  the  surface  a^ain,  tend  to  be  brought  together  so  that 
the  sound  energy  becomes  more  concentrated,  and  the  transmission 
loss  is  reduced  by  a  significent  amount.  Under  advantageous 
conditions,  these  rays  will  be  refocused  in  second,  third,  and 
higher  convergence  zones, 

(c)  Bottom  Bounce.  The  reflection  of  sound  off  the  bottom  and  back 
to  the  surface  provides  a  propagation  path  filling  in  the  ranges 
between  the  greatest  useful  range  of  the  mixed  layer  and  the 
innermost  edge  of  the  convergence  zone.  Unfortunately,  the 
unpredictable  and  usually  large  reflection  loss  that  occurs  at 
the  bottom  causes  the  transmission  loss  alon^  this  path  to  be 
rather  large. 
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ransmisrion  Loss  (dB) 


f  -  2  kHz 


a  *  0.08  dBAm 


e^ 


I  to  100  1,000 

Range  (km) 


(a)  Siirface  duct 

(b)  Convergence  zone 

( c )  Bottom  bounce 


Ref,  Urick,  Principles  of  Underwater  Sound  for  Engineers 
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CONVERCSI' ca  ZOI-i:  Pj\RAIIETERS 


|<-AC  -X 


Cq  •  the  maximuiii  speed  of  sound  above  the  SOFAR  channel  ajds 

A  c  »  the  difference  between  c^  and  the  speed  of  sound  at  the  elbow 

Z  »  the  vertical  distance  between  depths  for  which  c(z)  =  c^ 

r  _  =  the  ran£e  to  the  conver£ence  zone 

Snell's  Law  reveals  that  a  ray  which  is  horizontal  in  the  upper  ocean  VThere 
c(z)  ■  c  will  also  be  horizontal  at  a  distance  Z  deeper,  and  a  ray  which 
is  directed  downward  where  c(z)  ■  c^  must  penetrate  below  the  distance  Z. 

To  insure  that  enoufh  rays  are  returned  to  the  surface  to  focus  at  the 
convergence  zone,  the  total  depth  of  water  must  exceed  Z  or  Z  +  D  by  an 
amount  called  the  depth  excess.,  d 

If  there  is  no  surface  layer,  then  c  is  associated  vdth  the  siirface 
and  Z  is  measured  from  the  surface.  If  there  is  a  surface  layer,  then  c^ 
is  associated  with  the  bottom  of  the  layer  and  Z  measured  from  the 
bottom  of  the  layer. 


r*— AC 
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coi^TDKcacs  zom:  rai.c..';  m'd  irarsmij-'sion  loss 


A  rou£h  estimate  of  the  ran£e  at  vhich  a,  converfence  zone  forrns  is  £iven 

W  _ _ ^ 

J  _ 

2,8  Z  \Jc^/£^c  +  Layer  effect 


Layer  effect  »  0  if  there  is  no  surface  duct 

■  r  if  a  surface  duct  is  present  and  both  source 
a8d  receiver  are  near  the  top  of  the  layer. 


1'he  accuracy  of  the  above  expression  depends  on  hoi-j  veil  tue  constant- 
gradient  approximations  fit  the  true  c(z)  profile.  As  might  be  expected, 
vhen  the  eltox'J  of  the  curve  is  very  sharp,  and  the  gradients  above  and 
belov  the  elbox^  are  relatively  constant,  then  the  prediction  vill  be 
rather  good.  The  vidth  of  the  cor;verger’ce  zone,  over  vhich  the  rays 
are  appreciably  focused,  is  normally  about  10  per  cent  of  the  range  to 
the  zoi:e. 

Since  the  rays  begin  to  spread  out  spherically  as  they  leave  the 
source  and  are  attenuated  only  by  the  absorptive  losses  in  the  ocean, 
the  transmission  loss  for  convergence  zone  paths  can  be  written  in  the 
form 


TL  “  20  log  r  +  ar  -  G 


where  C  is  the  convergence  zone  gain.  Typical  values  of  G  range  between 
about  $  and  15  dB  when  a  convergence  zone  is  found. 
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Typical  Convergence  Zone  Kanges 


North  Pacific 


Surface  Temperature-”- 


Miniraum  depth  of  CZ  operation+ 


Range  to  first  CZ 


50°  F 

1270  ftm 

=  23fC  m 

i;7  kyd  = 

i 

1+3 

55 

1610 

29uO 

1 

'  52 

i+7 

60 

1900 

31+80 

<  56 

51 

65 

2150 

3930 

60 

55 

70 

2U00 

1+UOO 

61+ 

57 

75 

2600 

U760 

66 

60 

80 

2800 

5130 

69 

63 

!  Area 

1 

Surface  leraperature-it 

.'•linimum  depth  for  CZ 

1  operation+ 

1  1 

'  Ran(_e  to  first  CZ  ' 

1 

1 

1'!.  Pacific 

1  50°  F 

1  1270  ftm 

1  i 

f  U7  kyd 

1 

N.  Atlantic 

50 

i 

1  ^ 
i  1050 

1 

■  U6 

.  1 

Norwegian  Sea 

50 

1680 

53 

.•lediterranean  | 

67 

800 

33 

•it  Assumes  surface  duct  absent 
+  Allows  for  a  200  ftm  depth  excess 

Ref.  Convergence  Zone  Range  Slide  Rule,  Naval  tnderviater  Systems  ^enter 


44 


BOTTOM  HOUliCE 


For  an£les  of  depression  beyond  about  the  radius  of  curvature  of 

a  ray  is  so  lar^e  that  it  is  a  £ood  approximation  to  assume  Etrai£ht-line 
propagation  of  the  rays  over  the  distances  of  interest. 


A  ray  bouncing  off  the  bottom  ■will  return  to  the  surface  at  a  range 

r  =  2  n  tan  0 

Since  the  distance  travelled  by  the  ray  is  r/(cos  0),  and  since  the 
rays  spread  out  spherically  and  are  attenuated  only  by  the  absorptive 
losses  of  the  ocean  and  the  reflection  at  the  bottom,  an  approrimate 
transmission  loss  can  be  witten  as 


TL 


20  log 


r 

cos  0 


+ 


ar 

cos  0 


+  BL 


where  BL  is  the  bottom  loss  in  dBAounce. 


I 
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ECITOI'^  ]JDS£ 


Depan i? 


r/2H 


9 


Depends  on  bottom  composition 

Example:  At  21*  kHz  and  10°  £razin£  angle 


ilud 

Mud-sand 

£and-ittud 

Sand 

Stony 


16  dB/bounce 
10 
6 
h 
k 


Depends  on  surface  roughness 

Bottom  loss  goes  up  as  ratio  of  roughneSvS-to-viavelength  increases 


Ref.  Urick,  Pirinciples  of  Underwater  Sound  for  Engineers 
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Depth 


RSIIABLE  PROPAGATION  PATH  FOR  DEEP  SOURCE  OR  RECEIVER 


Reliable  Acoustical  Path  (RAP). 

If  the  target  is  shallow  and  the  sensor  deep,  or  vice  versa,  there 
exists  the  possibility  of  a  guarenteed  acoustical  path  between  the  two. 
This  path  is  useful  to  ranges  beyond  those  attainable  by  transmission  in 
the  surface  duct.  As  the  deeper  of  the  two  (sensor  or  target)  is  lowered 
further  and  further  below  the  layer,  the  shadow  zone  indicated  in  the 
diagram  moves  out  to  larger  ranges.  If  there  is  no  mixed  layer,  then  the 
leading  edge  of  the  shadow  zone  is  formed  by  the  downward-bending  ray 
which  just  grazes  the  surface. 

Because  the  angles  for  RAP  are  usually  fairly  large,  the  rays  can 
be  considered  to  travel  approximately  in  straight  lines  between  the 
target  and  sensor.  For  this  kind  of  propagation,  the  rays  spread  out 
spheiically  and  the  only  losses  are  those  from  the  absorption  of  sound 
in  sea  water.  An  approximate  transmission  loss  formula  is  then 


TL 


20  log 


r 

cos  6 


♦  ^  „ 
cos  0 


where  0  is  the  angle  of  elevation  of  the  rays  between  target  and  sensor. 

If  the  negative  gradient  (below  the  layer,  if  one  exists)  has 
magnitude  g,  then  the  range  to  the  shadow  zone  at  the  surface  is  estimated 
by  - 

l.li  VcQ(z-D)/g'  +  ^  Tg 

where  (z  -  D)  is  the  depth  below  the  bottom  of  the  layer.  If  there  is  no 
layer,  then  D  ■  0  and  there  is  no  skip  distance,  so  r^  ■  0. 
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Depth 


THE  SOFAR  CHAia'ffiL 


c(z) 


Range 


t 


V 


Just  as  for  the  surface  duct,  the  sound  spreads  out  more  or  less 
spherically  until  the  rays  are  bent  back  into  the  channel,  after  which 
the  sound  spreads  cylindrically.  Thus,  beyond  some  transition  ran£,e  r^* 
we  should  have  a  spreading  term  which  goes  as  10  log  r.  Those  rays 
whose  angles  of  inclination  with  respect  to  the  SOFAR  channel  axis  are 
large  enough  that  they  reflect  from  the  ocean  surface  or  bottom  will 
suffer  reflection  losses  and  scattering  losses  and  eventually  be  so 
weak  compared  to  the  trapped  rays  that  they  can  be  neglected.  Thus, 
for  propagation  in  the  channel,  there  are  no  losses  except  absorption, 
so  that  the  transmission  loss  has  the  simple  form 


TL  =  10  log  r  +  10  log  r^* 


for  ranges  beyond 


+  ar 


The  transition  range  depends  on  how  deep  the  sound  source  is  in  the  channel 
and  the  vertical  extent  of  the  channel.  For  sources  placed  reasonably 
down  in  the  channel,  the  transition  range  is  less  than  the  convergence 
zone  range  r  .  Notice  that  the  rays  cycle  up  and  down  in  the  channel, 
accomplisl'iing^one  cycle  over  an  interval  very  similar  to  the  convergence 
zone  range . 


A 
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TIME  STRETCIilr.G  IM  TlIC  SCF.'U'.  CiL\ra-,EL 


time 


L . 


An  interesting  fact  often  observed  in  the  SOFAE  channel  is  that  transient 
signals  are  "stretched"  in  time  as  they  propagate  to  large  ranges.  This  arises 
because  of  the  variation  in  the  speed  of  sound  over  the  vertical  extent  of  the 
channel.  For  the  SCF/VR  channels  commonly  encountered  in  the  northern  hemisphere, 
the  speed  of  sound  increases  sufficiently  rapidly  with  distance  above  or  below 
the  channel  axis  that  sound  traveling  oter  those  rays  which  swing  over  large 
vertical  distances,  like  a  and  b  in  the  above  sketch,  reaches  the  receiver 
at  times  earlier  than  sound  traveling  over  the  more  direct  paths  c,  d,  and  e. 
Thus,  the  sound  that  travels  straight  down  the  axis  of  the  channel  arrives  at 
the  receiver  last.  This  leads  to  the  characteristic  blossoming  of  the  received 
signal  from  an  explosive  charge,  the  sound  first  being'  very  weak,  and  then 
growing  in  volume  until  there  is  a  sudden  cutoff  of  sound  as  the  last  signals 
come  over  the  slowest  paths. 

The  exact  details  of  these  effects  depend  greatly '  on  the  character  of  the 
speed  of  sound  profile:  If  the  sound  does  not  increase  sufficiently  rapidly 
above  and  below  the  axis,  for  example,  the  time  stretching  may  not  occur,  or 
may  even  be  reversed,  vdth  the  sound  traveling  straight  down  the  axis  arriving 
first.  Additionally,  the  exact  location  of  the  source  and  receiver  with  respect 
to  tl'.e  channel  axis  is  very  important  with  respect  to  the  envelope  of  the 
received  signal  and  the  amount  of  time  stretching  encountered.  For  the 
tjipical  profiles  in  the  northern  henisphere,  the  increased  duration  i^t  of  the 
received  signal  over  that  sent  can  be  estimated  by 


At  .  1  Ac 

t  3  c 


if  source  and  receiver  are  on  the  channel  axis,  and  less  if  they  are  not. 
The  quantity  t  is  the  total  time  of  flight  between  source  and  receiver  and 
Ac  is  the  maximiuti  variation  in  the  speed  of  sound  between  the  upper  and 
lower  boundaries  of  the  channel  indicated  in  the  previous  page.  Typically, 
this  leads  to  about  9  seconds  of  stretching  for  each  1000  nautical  miles 
the  signal  has  travelled. 
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NORMAL  MODES 


Another  approach  to  the  propagation  of  sound  in  the  ocean  is 
Normal  Mode  Theory.  If  the  mathematical  difficulties  of  this 
approach  can  be  overcome,  it  yields  more  exact  values  for  the 
transmission  loss  than  can  be  obtained  using  Ray  Theory. 

Normal  Mode  Theory  yields  a  sum  of  individual  solutions 
(each  called  a  normal  mode)  which  must  be  combined  to  obtain 
a  complete  expression  for  the  sound  field. 


Each  one  of  these  normal  modes  has  its  own  behavior 


This  part  looks  just  like  a  cylindrical  wave 
traveling  outward  with  its  amplitude  falling 
off  like  l/Zr' .  Notice  however  that  it  has 
its  own  special  speed  c^. 


This  pairt  is  an  amplitude  factor  which  changes  with 
depth.  Notice  that  each  normal  mode  has  its  own 
special  depth-dependent  amplitude  factor. 


These  special  factors  are  called  "normalized  eigenfunctions" 
They  are  found  by  solving  the  wave  equation  and  specifying  the 
behavior  of  the  pressxire  at  the  ocean  surface  and  bottom.  The 
solution  for  each  also  gives  a  required  value  for  c^.  Both 
the  depth  dependence  of  Z^  and  the  value  of  Cj^  depend  on  the 
frequency  of  the  soxxnd  waves. 

This  kind  of  a  solution  can  be  obtained  whenever  the  speed  of 
sound  profile  c(z)  is  shaped  so  that  the  sound  can  be  trapped 
in  a  channel.  The  trapping  of  sovind  in  the  mixed  layer  is  one 
such  example,  and  the  propagation  of  sound  to  long  ranges  in 
the  SOFAR  channel  is  another. 
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NORMAL  MODES  FOR  A  CHANNEL 


The  above  sketches  represent  the  normalized  eigenfvinctions  of  the 
first  three  and  sixth  normal  modes  for  so\xnd  trapped  in  the 
SOFAR  channel.  These  sketches  show  that  the  greatest  part  of 
the  pressure  field  is  found  in  the  channel,  but  the  long  "tails" 
which  extent  to  depths  far  below  the  channel  axis  show  that 
there  is  sound  ensonification  at  depths  greater  than  those 
predicted  by  drawing  ray  paths. 
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CUTOFF  FREQUEiNCY  AND  LOSSES 


1 

» 


1 


Each  normal  mode  has  its  own  cutoff  frequency  f  .  If  the  sound 
has  a  frequency  below  that  of  cutoff  for  a  particular  mode,  then 
that  mode  cannot  be  excited  to  carry  so\ind  energy  in  the  channel. 

If  the  frequency  of  the  soiond  is  slowly  increased  with  time,  the 
above  sketches  show  how  a  normal  mode  can  be  excited.  Just  above 
cutoff,  the  normal  mode  has  a  very  long  tail  which  goes  to  great 
depths  beyond  the  edges  of  the  channel.  As  the  frequency  increases 
further  above  cutoff,  the  tail  gets  shorter  and  shorter,  until 
at  frequencies  well  above  f  there  is  very  little  soiind  existing  • 
outside  the  channel. 

This  means  that  for  a  specified  frequency,  souhd  can  be  trapped  and 
sent  to  large  distances  only  in  those  normal  modes  whose  cutoff 
frequencies  are  less  than  the  frequency  of  the  sound.  Normal  modes 
with  higher  cutoff  frequencies  will  have  longer  tails  than  those 
with  lower  cutoff  frequencies,  and  those  whose  cutoff  frequencies 
exceed  the  frequency  of  the  sound  cannot  be  excited  at  all. 

The  larger  the  tail  extending  out  of  the  channel,  the  more  energy 
can  be  lost  from  the  sound  in  the  channel.  This  means  that 
higher  modes  (those  with  higher  cutoff  frequencies)  will  attenuate 
faster  and  be  less  significant  at  great  distances.  This  effect  is 
sometimes  called  mode  stripping. 

These  energy  losses  for  each  mode  can  be  handled  by  specifying  an 
absorption  coefficient  for  each  normal  mode, 

where  increases  with  increasing  n. 


Then,  we  have 

e 
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TRANSMISSION  LOSS  (NORMAL  MODES) 


/L 


If  incoherence 

P=\/p.’ 


low  frequency 
short  ranges 
homogeneous  water 
smooth  boundaries 


high  frequency 
long  ranges 
inhomogeneous  water 
rough  boundaries 
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Transmission  Loss  in  dB 


RAYS ,  NORMAL  MODES ,  AND  REALETY 

Normal  mode  theory  is  a  complicated  mathematical  formalism  which  does  not 
lend  itself  to  simple  intuitive  interpretation.  However,  it  is  an  exact 
theory  which  avoids  certain  approximations  inherent  in  ray  theory.  It  is 
useful  in  discribing  situations  for  \diich  ray  theory  fails.  An  example 
of  one  such  situation  is  illustrated  below. 


Range 


The  soxirce  is  in  the  layer  and  the  receiver  is  below  the  layer,  Ray 
theory  predicts  a  shadow  zone  at  large  enough  ranges,  VJhile  there  can 
be  some  scattering  from  the  ocean  surface  into  the  shadow  zone,  this 
will  be  small  for  low  frequencies  and  low  sea  states.  The  surface 
interference  pattern  at  short  ranges  shows  that  the  sea  was  quite 
smooth  (the  observed  sea  state  was  about  1  or  2). 


Notice  that  the  experimental  data  show  significant  ensonification  of  the 
shadow  zone.  This  cannot  be  predicted  by  the  ray  theory,  and  the  scattering 
of  sound  from  the  surface  is  not  strong  enough  to  account  for  it. 

Reference:  Pederson  and  Gordon,  Jour,  Acoust.  Soc.  Am.  105  (1965) 
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SWL4IIY  OF  TRANSMISSION  LOSS  EXAITLES 


For  both  source  and  receiver  in  the  surface  layer,  and  at  ranj^es 
£reater  than  the  transition  range,  the  transmission  loss  can  be 
estimated  by 


TL  “  10  log  r.  +  10  log  r  +  ar  +  br/r  . 

t  s 


For  a  convergence  zone  detection,  we  can  v’rite 
TL  “  20  log  r  +  ar  -  G  , 


If  the  propagation  path  is  bottom  bounce,  and  the  angle  of  depression 
of  the  sound  beam  in  that  mode  is  sufficiently  large,  then  the  simple 
geometrical  model  gives 


TL 


20  log 


r 

cos  0 


+ 


ar 

cos  9 


+ 


BL  . 


For  propagation  betvjeen  a  deep  point  and  a  shallow  point  over  the 
reliable  acoustical  path,  an  approximate  transmission  loss  formula  is 


TL 


20  log 


ar 


cos  9 


cos  9 


Finally,  for  the  propagation  of  sound  trapped  in  the  SOFAR  channel, 
if  the  range  is  greater  than  the  transition  range  tbe  transmission 
loss  can  be  estiraated  by 


TL 


10  log  r  +  10  log  r^’ 


+  ar  . 
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COMBININS  SIGNALS 


Now,  let  us  no  longer  restrict  ourselves  to  tonals  {monofrequency 
signals).  Because  we  really  don't  deal  with  only  a  single  frequency, 
we  have  to  find  out  how  signals  of  different  frequencies  combine. 
Consider  the  combination  of  a  large  number  of  monofrequency  signals* 


The  combination  of  many  signals,  each  a  pure  tonal,  will  result  in 
a  total  signal  which  does  not  appear  to  have  definite  amplitude 
or  period.  However,  the  composite  signal  is  an  acoustical  wave, 
and  the  total  acoustical  energy  it  possesses  is  the  sum  of  the 
acoustical  energy  of  each  of  the  individual  tonals.  Thus, 
if  E<r  is  the  total  energy  of  the  composite  signal,  we  have 

N 

Et  =  El  +  £3  +  •  •  •  =  Ei 

i=  I 

where  E^  is  the  energy  of  signal  1,  £3  that  of  signal  2,  and  so  on. 

For  the  complicated  signals  encountered  in  real  life,  which  have 
large  variations  from  moment  to  moment  in  their  amplitude  and 
shape,  the  frequencies  of  the  individual  tonals  which  combine  to 
form  them  are  extremely  closely  spaced.  Indeed,  it  is  often 
more  convenient  to  talk  about  the  energy  contained  in  the 
collection  of  tonals  which  lie  between  some  specific  frequency  and 
another  frequency  which  is  only  inf initisimally  different. 

If  aE  is  the  total  energy  of  all  the  tonals  lying  within  the 
incremental  bandwidth  aw,  then  we  can  define  the  quantity  £ 
as  the  average  energy  density  over  the  incremental  bandwidth  aw  by 


A.MJ- 


Notice  that  the  energy  density  has  the  units  of  energy/Hz.  If  we 
analyze  some  complicated  signal,  therefore,  we  can  describe  it 
in  terms  of  the  energy  density  it  contains  in  each  incremental 
interval  of  frequency  over  the  entire  range  of  frequencies  which 
comprise  it. 


N 


•  > 

'\AAAA.' WyW>J 
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ENERGY  SPECTRUM  AND  PRESSURE  SPECTRUM  LEVEL 

An  example  of  the  analysis  of  some  complex  signal  into  the 
energy  densities  found  in  each  frequency  interval  is  suggested 
in  the  following  curve i 


The  total  energy  of  the  signal  is  given  by 

I 

which  is  the  area  under  the  curve  of  energy  density  vs  frequency. 

If  all  of  the  are  chosen  to  be  1  Hz  wide,  then  the  graph  gives 
the  energy  to  be  found  in  each  1  Hz  interval,  and  the  curve  is 
called  the  energy  spectrum  of  the  signal. 

It  is  more  conventional  in  acoustics  to  describe  the  signal  not 
in  terms  of  the  energy  spectrum,  but  rather  in  terms  of  the 
pressure  spectrum  level  PSL,  which  we  shall  define  here  as 

VSL-.  2  0  JBrs 

where  the  quantity  p  ,  the  effective  pressure  per  Hz,  generates 
the  associated  energy  found  in  the  1  Hz  bandwidth.  For  example, 
if  we  receive  a  certain  acoustical  signal,  pass  it  through  a 
filter  with  a  1  Hz  bandwidth,  and  obtain  the  effective  pressure 
amplitude  of  the  resultant  filtered  signal,  say  100  ^Pa,  then  the 
value  of  p  for  that  1  Hz  bandwidth  is  100  /ixPa/l  Hz  *  100  ^Pa/Hz. 

A  representative  curve  of  the  pressure  spectrum  level  vs  frequency 
is  presented  below t 
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AVERAGED  ENERGY  SPECTRUM  AND  AVERAGED  PSL 


Very  often,  what  is  presented  as  the  energy  spectrum  or  the 
pressure  spectrum  level  is  not  that  which  would  be  measured  with 
a  1  Hz  filter,  but  is  rather  calculated  from  measurements  made 
with  a  filter  of  larger  bandwidth. 

The  overall  energy  to  be  found  in  some  arbitrary  bandwidth 
w  is  the  area  under  the  energy  spectrum  vs  frequency  curve 
over  the  appropriate  range  of  frequencies; 


This  gives  the  total  energy  found  within  this  bandwidth,  and 
does  not  say  how  e  is  distributed  over  w.  We  can  therefore  write 
the  averaged  energy  spectrum  £  as  £  =  E^/w.  A  plot  of  £  vs 
frequency  would  give  a  smoothed  curve,  lacking  the. detail  of  the 
true  curve,  but  sufficient  to  calculate  the  total  energy  present 
over  bandwidths  comparable  to  or  larger  than  the  bandwidths  selected 
to  obtain  g"  . 

Analogously,  it  often  happens  that  what  is  presented  as  a 
pressure  spectrum  level  is  in  reality  a  smoothed  curve,  the  result 
of  plotting  data  obtained  with  filters  having  bandwidths  appreciably 
larger  than  1  Hz  t 


Use  of  too  broad  a  bandwidth  in  calculating  averaged  spectrum 
levels  will  wash  out  details.  The  spikes  in  the  true  spectrum 
level  in  the  above  example  are  lost  to  view  because  the  bandwidth 
over  which  the  averaged  levels  are  taken  is  so  wide  that  the  energy 
from  the  spike  is  small  compared  to  all  the  energy  from  the 
background  in  each  bandwidth. 
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THREE  WAYS  OF  DOING  IT 


V 


I. - A— 

Frequency 


The  first  graph  presents  the  true 
spectrum  level  and  the  other  three 
present  the  same  spectrum,  except 
that  the  spectrum  levels  have 

been  averaged  over  the  bandwidth 
w.  Graph  2  presents  the 
averaged  spectrum  level  in  a 
bar  graph  format.  Graphs 
3  and  4  connect  the  data 
with  either  straight-line 
segments  or  a  smooth 
.4>  curve.  Notice  that 
Graph  4  gives  no 
information  relating  to 
the  bandwidth  of  the  filter 
used  to  obtain  the  averaged  spectrum 
levels?  this  would  have  to  be 

indicated  separately,  and  is  not 
likely  to  be  done.  Graph  3  is 
somewhat  better,  but  the  use 
of  the  straight-line  segments 
seduces  the  eye  into  assuming 
that  the  true  spectrum  level 
j  closely  follows  the 

j  straight-line  segments. 

I  Graph  2  is  the  least 

subject  to  being 
misinterpreted  I  The 
width  of  the  individual 
bars  gives  the  bandwidth 
of  the  filter  quite 


Frequency 


visually,  and  each  bar  strongly 

emphasises  by  its  flat  top  that  the 
value  of  the  PSL  appropriate  to 
that  frequeYicy  interval  is  the 
averaged  value. 

Because  of  the  fact  that 
the  averaged  spectrum  level 
smooths  out  and  erases 
detail  which  would 
otherwise  be  found  in 
the  true  spectrum  level, 
it  is  very  important 
that  the  extent  of  averaging 
over  frequency  be  well- 


Frequency 


indicated.  This  is  accomplished  quite 
easily  and  unambiguously  by  the  format 
of  either  Graph  2  or  Graph  3*  Further, 
Graph  2  clearly  indicates  the 
averaging  of  data  by  the  flat-topped 
nature  of  the  curve,  whereas 
Graph  3  can  tend  to  mislead 

the  eye  in  a  casual  examination. 


Frequency 


SPECTRUM  LEVEL  AND  BAND  LEVEL 


Recall  that  the  total  energy  in  a  bandwidth  w  was  the  area 
under  the  energy  density  vs  frequency  curve.  The  analogous 
quantity  in  terms  of  pressure  levels  is  the  band  level  BL.  This 
is  the  expression  in  dB  terms  of  the  total  effective  pressure 
generating  the  acoustical  energy  in  the  specific  bandwidth. 

NOTICE  THAT  THE  BAND  LEVEL  IS  NOT  THE  AREA  UNDER  THE 
SPL  vs  FREQUENCY  CURVE. 


This  is  because  the  log  of  a  product  is  the  sum  of  the  logs  of 
the  individual  factors. 

Consider  a  plot  of  energy  density  £  vs  frequency  as  showni 


£  > 

t 

o 

k 

AaT  - 

_ 

f 

It  is  clear  that  the  total  energy  E  in  the  bandwidth  w  is  given  by 

E  =  £„w 

If  we  express  this  in  log  form,  we  have 


It  turns  out  in  acoustics  that  taking  10*  log(E/E^gj.)  is  the  same 
as  taking  20*  log(P/Pj.gf ) .  This  means  that  the  above  equation 

is  the  same  as 


ZoM 


-  20  \oJl. 


1  N' 


or,  treating  the  quantity  1  Hz  as  understood  but  suppressed. 


BL  =  PSL  +10  log  w 

In  general  this  expression  is  exact  only  if  the  spectrum 
level  is  constant  over  the  bandwidth  w.  However,  it  is 
reasonably  accurate  if  the  true  spectrum  level  does  not 
fluctuate  by  more  than  a  couple  of  dB  over  the  range  of  w. 


Notice  that  this  formula  can  be  turned  around i  If  a  band 
level  BL  has  been  experimentally  obtained  by  measuring  the 
acoustical  energy  over  some  bandwidth  w,  then  the  averaged 
spectrum  level  over  this  bandwidth  can  be  found  from 

PSL(averaged )  =  BL  -  10  log  w 
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COMPARISON  OF  A  SPECTRUM  LEVEL  AND  A  BAND  LEVEL  CURVE 


/  V 
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Because  of  the  (10  log  w)  term  which  must  be  added  to  the 
PSL  to  obtain  the  appropriate  BL  for  a  certain  frequency  interval, 
it  is  seen  that  the  shape  of  the  BL  curve  is  such  that  the  BL 
is  stronger  at  higher  frequencies.  This  is  because  the  BL  measures 
the  total  acoustical  energy  in  some  frequency  interval,  whereas  the 
PSL  measures  the  energy  per  unit  Hz. 

Both  kinds  of  curves  have  their  advantages  and  disadvantages. 

A  curve  of  the  true  spectrum  level  gives  the  maximum  possible 
amount  of  information,  but  requires  considerable  mathematical 
analysis  if  it  is  desired  to  obtain  the  band  level  within  some 
specified  frequency  interval  from  the  spectrum  level.  A  curve 
of  averaged  spectrum  level  is  not  quite  so  difficult  (timeconsuming) 
to  use,  and  band  levels  from  it  can  be  found  more  easily,  but  often 
important  details  are  lost.  A  curve  of  band  level  is  quite  direct 
to  use,  as  long  as  the  bandwidth  agrees  with  the  particular  receiving 
system  being  used,  but  it  also  has  smoothed  out  much  detail. 
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A  SAMPLE  CALCUIATION 


As  a  simple  example,  let  us  obtain  the  overall  band  level 
for  the  PSL  curve  presented  below t 


PSL 

•Q  f 


f 


From  the  steplike  nature  of  the  curve,  we  can  specify  two  bandwidths 
w^  and  W2  as  shown,  with  the  associated  averaged  spectrum  levels. 

For  each  of  these,  we  can  calculate  the  appropriate  band  level i 


BL(wj^ )  =  B(w^)  +  10  log  Wj^ 
BL(w2)  =  B(w2)  +  10  log  W2 


So  far,  we  have  had  to  do  no  more  than  a  little  curve-fitting  and 
application  of  formulas  which  have  been  presented  previously. 

Now,  however,  a  new  problem  has  appeared.  We  must  combine  the 
two  band  levels  into  a  single  bandlevel  for  the  total  frequency 
interval  w^^  +  ^2*  combined  BL  is  not  the  sum  of  the  individual 

ones,  because  we  are  in  dBi  Energies  add  but  logs  of  energies 
do  not.  What  must  be  done  is  to  determine  the  amount  of  energy 
given  by  BL(wj^ )  and  the  amount  of  energy  given  by  BL(w2),  and 
these  energies  added  together  and  then  expressed  as  a  band  level. 

Let  us  illustrate  this  with  an  example.  We  shall  first  do 
the  example  in  the  method  indicated  above,  and  then  present  a 
simple  nomogram  which  avoids  all  the  work. 


13  ('*^1  )-  o  i  ^ 
=  I  35" 


x-r,  r  1000 
‘~LOOO 


BL(. 


10'^'^)  -  112.1  ^ 
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NOMOGRAM  FOR  COMBINING  LEVELS 


As  an  example  of  the  application  of  this  nomogram,  recall  the 
problem  worked  out  on  the  previous  page.  The  two  band  levels 
which  were  calculated  there  were  BLj^  =  170  dB  re  l/^Pa  and 
BLo  =  168  dB  re  IjLAPa.  The  difference  is  +2  dB;  finding  this 
value,  projecting  upward,  and  reading  to  the  left  gives  2.1  dB, 
so  that  the  combination  of  these  two  levels  yields  a  total 
band  level  of  170  +  2.1  =  172.1  dB  re  l^Pa.  Since  the  usual 
uncertainty  in  levels  exceeds  a  few  tenths  of  a  dB,  it  is 
appropriate  to  round  this  combined  level  to  172-  dB  re  l^xpa. 


TONALS  COMBINED  WITH  BACKGROUND 


-J 


-.'jf  >■ 


Frequency 


When  there  are  strong  tonals  present,  the  band  level  cannot 
be  expressed  as  B  +  10  log  w,  because  of  the  additional  energy 
contributed  by  the  spike(s).  We  must  calculate  the  energy  from 
the  background  spectrum  (without  the  spikes)  and  add  this  to  the 
energy  of  the  monofrequency  spikes  which  are  also  present: 

BL3  »  B  +  10  log(w-l )  =  B  +  10  log  w  if  w  ^  1 

BL^  =  A  +  10  log  1  =  A. 

The  total  band  level  is  found  by  combining  BL^  and  BLg  with  the 
help  of  the  nomogram. 

Notice  that  if  the  tonal  is  relatively  strong  and  the  bandwidth 
narrow  enough,  the  combined  level  will  be  due  essentially  to  the 
tonal.  As  the  bandwidth  w  increases,  however,  the  influence  of 
the  (10  log  w)  term  will  become  stronger,  causing  BLg  to  increase. 
Eventually  the  contribution  from  the  background  will  get  sufficiently 
large  that  BLg  will  be  sufficiently  large  that  the  influence  of 
BL^  will  be  lost. 
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TIK  SONATi 


I  EL  ;:iML  +  DT 

I 


EL  =  Echo  Level  for  active  sonar  or  the  received  Signal  Level  for  passive 
sonar.  It  is  the  level  in  dB  of  the  desired  si£nal  vhen  it  reaches 
the  receiver. 

IIL  =  Masking  Level.  It  is  the  level  in  dB  of  the  undesired  signal  at  the 
receiver  v?hich  is  competing  vdth  the  Echo  Level. 

DT  *  Detection  Threshold.  This  is  the  number  of  dE  by  vhich  the  Echo  Level 
must  exceed  the  Masking  Level  to  allov;  the  receiving  system  to 
register  a  detection  ■with  a  specified  degree  of  confidence.  This 
quantity  involves  the  probability  of  detecting  the  target  in  the 
presence  of  noise  and  also  the  probability  of  gi'ving  a  false  alarm. 


The  SCliAE  equation  states  that  signal  at  the  receiver  must  equal  or  exceed 
background  noise  by  an  amount  specified  by  the  Detection  Threshold  if  a 
specified  degree  of  performance  is  to  be  achieved. 
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LIST  CF  Sr.-BOIB 


EL 

m 

Echo  Level 

ML 

m 

Masking  Level 

DT 

m 

Detection  Threshold 

SL 

m 

Source  Level 

TL 

m 

Transmission  Loss 

PSL 

m 

Pressure  Spectrum  Level 

BL 

m 

Band  Level 

SSL 

m 

Source  Spectrtim  Level 

AG 

m 

Array  Gain 

NL 

m 

Noise  Level 

NSL 

m 

Noise  Spectrum  Level 

DI 

m 

Directivity  Index 

w  *  bandwidth 

d  -  detection  index 

wt  «  time-bandwidth  product 
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PASfUvii  soiim 


o 

EL  =  SL  -  TL 


SL  -  TL  ^  ML  +  DT 


vhere  SL  »  Source  Level  of  the  target  . 
TL  =  One-way  transmission  loss 


Each  term  in  the  above  equation  vdll  nov?  be  discussed  separately. 
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(s^  -  TL  ^  ML  +  DT 
SOURCE  LEVEL  -  Pj\S?IVE  SOLAR 


.  Tonals  (monfrecuency  coraponants) 

Sources:  Propulsion  machinery  (Main  motors,  reduction  tear) 

Auxiliary  machinery  (generators,  pumps,  air-conditioners) 
Propeller-induced  resonant  hull  excitation 
I iydrodynamically -induced  resonant  excitation  of  cavities, 
plates,  and  appendaj^es 


SL  =  Level  of  acoustic  energy  1  meter  from  tie  acoustic  center  of  the 
source. 

i’Jeasurements  are  usually  made  at  a  distance  from  tlie  source  and 
extrapolated  back  to  1  meter. 


II .  Hroadband  noise 

Sources:  Cavitation  at  or  near  the  propeller 
Raiiiated  flow  noise 
Caviation  at  struts  and  appendages 


SL  »  SSL  +  10  log  w 

where  w  is  the  bandid-dth  of  the  system  and 

SSL  =  pressure  spectrum  level  1  meter  from  the  acoustic  center 
of  the  source. 

If  the  SSL  is  not  essentially  constant  over  w,  the  bandwidth 
must  be  subdivided  into  smaller  segments,  the  LL  of  each 
found,  and  then  all  the  BL's  combined  to  give  SL 
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TIPI CAL  SOURCE  SPECTRUM  LEVELS 


lypical  source  spectrum  levels  for  a  conventional  submarine  on  electric 
drive  at  periscope  depth.  SSL  re  l^Pa. 


Adapted  from  Urick,  Principles  of  Underwater  Sound  for  Engineers 


SL  -  TL  nr 

INKING  LEVEL 


Unwanted  noise  may  come  in  uniformily  from  all  direction  or  it  may  arrive 
from  a  pr^fered  direction. 


ML  -  NL  -  AG 

where 

ML  ■  raaskin£  ictel 

NL  ■  noise  level  in  the  environment 

AG  “  array  gain,  the  improvement  (reduction)  in  dB  against  the 
noise  field  resulting  from  the  array  design. 
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where 


and 


Kote: 


I4L  «  NL^  -  AG 
Koise  Level 


NL  =  NSL  +  10  lot  V 


NL  »  noise  level,  the  acoustic  energy  in  dB  found  in  the  bandwidth 

of  the  system  ^ 

NSL  »  noise  spectrum  level,  pressure  spectrum  level  of  the  unwanted 

noi se , 


w  *  the  bandwidth  of  the  system  in  ilz. 


The  same  limitations  apply  in  using  this  formula  as  were 
pointed  out  in  the  discussion  of  the  broadband  SL, 
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NL  *  :  WSIJ  +  10  log  w 


NOISE  SPECTRUM  LEVEL  (Deep  water) 


Adapted  from:  Venz,  Jour.  Acoust.  Soc,  Am,  3U  ,  1536-1956  (1962) 
Perrone,  ibid  i+6,  762-770  (196^ 
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rJL  =  (USL  +  10  l0[  w 
i  SP-iCTKlM  LCViL  (Self -noise) 


Sources  of  self-noise:  Hydrodynamic  noise  due  to  flov  past  receiver 

Acceleration  of  receiver 
Machinery  noise  of  ov;n  "ship" 


Values  usually  presented  as  ai;  "equivalent  omnidirectional"  noise  spectrum 
level. 


KL(  self -noise )  »  hSLCself-i'Oise)  +  10  lot  v 


If  sipnificant  self -noise  is  present,  it  must  le  combined  vdth  the  ambient 
noise  level  to  obtain  the  total  noise  level.  Since  these  two  sources  of 
noise  are  independent,  their  snerpies  add  and  the  nomopram  can  be  used. 
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ML 


NL  -fAG 


ARRAY  GAIN  AND  DIRECTIVITY  INDEX 


Array  gain  is  the  decrease  in  masking  level  from  the  noise 
level  brought  about  by  the  design  of  the  receiving  array. 


If  the  noise  is  isotropic  (arriving  uniformly  from  all  directions) 
and  the  signal  arrives  from  one  direction  only,  then 

AG  =  DI 


where 


DI  =  Directivity  Index 

=  10  log  (the  ratio  of  the  noise  energy  sensed 

by  an  omnidirectional  receiver  to  that 
sensed  by  the  directional  receiver) 


In  this  case, 


ML  =  NL  “  DI 


Array  gain  can  include  both  the  design  of  the  receiving  array 
and  also  the  signal  processing.  Signal  processing  such  as  correlation, 
by  taking  into  account  the  different  spatial  or  temporal  properties 
of  the  signal  and  noise  can  result  in  an  array  gain  different  from 
the  directivity  index. 
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Pair  hj  pair, 
the  signals 
cancel 


axial ^ 


DI  -  0  dB 


Onnidirectional  receirer 


m  •  20  dB  Receiflng  only  1$  of  the  energy  present  in  the  noise 
field. 

I  depends  on  sise  of  the  array  compared  to  the  vave length.  For  example , 
or  a  long  line  array  of  length  L  such  greater  than  ^ 

nr  -  10  log  (2L/A  ) 

snerally,  as  frequency  increases  the  receiver  becomes  more  directive  (m  Increases) 
lit,  if  the  frequency  is  too  high  we  may  get  more  major  lobes. 
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A  PHYSICAL  MEANING  OF  DIRECTIVITY  INDEX 


Surface  area  of  the  unit  sphere  »  Htt  or 

Area  of  the  illuminated  patch  ■  _Os.O)' 

where  JX  is  the  solid  angle  subtended  by  the  patch. 

The  fractional  area  of  the  sphere  which  is  illuminated 
is  given  by  JX  /4'Tr  and,  for  this  simple  example,  the  directivity 
D  is  given  by  ^ 

The  directivity  index  is  the  expression  of  the  directivity  D  in 
dB  terms, 


Di  =■  /o  D 


A  more  general  definition  of  the  directivity  is  given  as 
follows:  If  the  pressure  field  generated  by  an  array  acting  like 

a  source  (or  the  relative  sensitivity  of  an  array  acting  like  a 
receiver)  is  written  in  the  form 

P(r,0,^)  -  Pax(r)  ii(9,0) 

where  P^x  is  the  pressure  on  the  acoustical  axis  of  the  array, 
then  H(0,j^)  is  normalized  to  have  a  maximum  value  of  unity  on 
the  acoustical  axis,  and  to  be  no  greater  than  unity  for  any  other 
direction.  The  directivity  can  then  be  expressed  as 
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SOME  ROUGH  ESTIMATIONS  OF  DI 


If  we  have  a  line  array  with  a  single  major  lobe,  and  if 
the  angular  width  of  the  major  lobe  is  9  (in  radians),  then  if 
the  array  is  surrounded  by  a  unit  sphere,  the  illuminated  patch 
appears  as  below 


and  the  directivity  can  be  crudely  estimated  to  be 

- 2. 


2n(iW) 


e 


(  9  in  radians  ) 


where  it  must  be  understood  that  9  is  small. 


A  little  more  careful  estimation  of  the  directivity  from 
the  geometry  of  the  radiation  (or  reception)  pattern  leads  to 
some  relatively  simple  but  only  approximate  estimates  for  the 
directivities  of  arrays  of  simple  design. 

Line- like  arrays: 


r>  ^  I 

^  I  A 


»  I 


Piston-like  arrays: 


»l 


These  formulas  can  be  used  as  rough  estimates  if  the  arrays  are 
highly  directive  (DI  in  excess  of  about  6  dB  at  the  least),  are 
designed  to  have  only  one  major  lobe,  and  are  not  strongly  shaded 
or  steered. 

An  array  is  shaded  when  different  elements  have  different  amplifications. 
This  is  done  to  alter  the  width  of  the  major  lobe  and  the  relative  strengths 
of  the  side  lobes.  Narrowing  the  major  lobe  will  increase  the  strengths  of 
the  side  lobes,  and  vice  versa.  An  array  is  steered  to  point  the  major  lobe 
in  various  directions  with  respect  to  the  array.  This  is  done  by  introducing 
time  delays  in  the  outputs  of  the  elements. 


77 


PASSIVE  SONAR  EQUATIONS 


Let's  put  it  all  together: 

BROADBAND  WITH  NEGLIGIBLE  TONAIS 

SSL  +  10  log  V  -  TL  >  NSL  *  10  log  w  -  AG  ♦  DT 

SSL  -  TL  ^  NSL  -  AG  ♦  DT 


TONALS  WITH  NEGHGIBLE  CONTRIBUTION  FROM  BACKGROUND 


SL  -  TL  >  NSL  +  10  log  w  -  AG  ♦  DT 
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FEGURE  OF  MERIT 


(POM) 


The  passive  SONAR  equations  can  be  rewritten  by  solving  forTL: 


TL  ^  POM 

where 

POM  -  SL  -  (NSL  +  10  log  w  -  AG  +  DT) 
or 

POM  -  SSL  -  (NSL  -  AG  +  DT) 
for  tohal  or  broadband  detection  respectively. 
Graphically: 


ranges  for  which  detections  are  possible  at  conditions 
equal  to  or  better^  than  specified  by  DT, 
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To  take  advantage  of  a  narrow  bandwidth  and  still  permit  reception  of  a 
Doppler  shifted  signal,  it  is  necessar7  to  design  a  system  that  en^iloys 
many  parallel  filters  irtiich,  while  individtially  narrow,  cover  the  entire 
frequency  range  of  interest. 

The  output  of  ouch  a  system  might  be  presented  in  this  manner: 


At  any  instant,  the  output  of  each  of  the  individual  filters  might  look 
like  this: 
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Probability  Density  of  An5)litude  A 


SL  -  TL  <  ML  +  @ 
DETECTION  THRESHOLD 


The  detection  process  consists  of  designating  a  threshold  which,  when  exceed, 
causes  a  detection  to  be  recorded*  If  the  signal  is  much  stronger  than  the 
noise,  it  is  clear  that  a  threshold  can  be  defined  that  will  allow  valid 
signals  to  be  recorded  while  '  ignoring  the  noise.  However,  when  the 

signal  and  noise  are  of  comparable  size,  any  threshold  that  will  catch  a 
reasonable  number  of  valid  signals  will  also  record  "detections"  when  a  valid 
signal  is  absent* 

The  probablity  of  obtaining  a  certain  an^^litude  of  the  noise  alone  as  well  as 
for  the  signal- and' noise  is  sketched  below: 


where 

Aj^  is  the  mean  value  of  the  amplitude  with  noise  alone 
A  is  the  mean  value  of  the  amplitude  with  signal  and  noise 
<r  is  a  measure  of  the  width  of  the  curves  (assumed  equal) 

T  is  the  threshold 

The  shaded  area  marked  "detection"  is  the  probablity  of  making  a  valid  detection 
when  the  threshold  is  T.  (The  area  under  the  curve  is  taken  as  unity*)  And 
the  shaded  area  marked  "false  alarm"  is  the  probablity  of  making  a  "detection 
when  the  signal  is  absent* 

Let  P(D)  •  the  probability  of  detection 

P(FA)  ■  the  probability  of  a  false  alarm 
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DETECTION  INDEX 


An^jlitude 


We  need  to  find  some  simple 
criterion  which  will  describe 
how  the  cixrves  of  noise  alone  and 
signal  and  noise  combined  overlap, 
since  for  a  given  value  of 
threshold,  it  is  the  overlap  of 
the  curves  which  establishes 
the  relative  values  of  P(D)  and 
P(FA). 

Notice  the  configuration  of  the 
two  curves  given  in  the  first 
graph.  The  curves  are  quite  well- 
sepsirated,  in  that  their  peaks 
have  amplitudes  idiose  difference 
is  larger  than  the  o'  .  These 
two  curves  could  be  considered 
well-eeparated . 

In  the  second  graph,  each  curve 
has  the  same  average  value  of 
A  as  before,  but  now  the  0*^13 
twice  as  large.  The  curves  now 
are  not  well-separated. 

In  the  third  graph,  the  curves 
have  the  same  O''  as  in  the  second, 
but  the  two  have  been  shifted 
further  apart  so  that  they  are 
well-separated.  Indeed,  con5)ari3o: 
of  the  first  and  third  graphs 
will  reveal  that  the  amount  of 
overlap  of  the  pairs  of  curves 
is  the  same.  Notice  that  the 
ratio  of  the  difference  in  the 
verage  values  of  the  aunplitudes 


— f - 1 - - i- 

A  y  A 


^s&n  ^  ^  ^ 

first  graph  is  identical  with 

that  for  the  third.  Thus,  it 

is  the  separation  between  the 

two  curves  ih  units  of  the  width 

of  the  curves  which  measures 

how  well-separated  the  two  curves 

are.  It  is  thus  plausible  to 

define  the  detection  index  d  by 
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Let  us  study  the  effects  of  changing  the  positions  of  the  threshold  and  the 
detection  index  on  the  probabilities  of  detection  and  faCLse  alarm.  If  the 
original  configuration  of  the  rusise  curve,  signal  and  noise  ciirve,  and 
threshold  are  as  given  in  the  top  curve,  then  we  can  establish  the  indicated 
P(D)  and  P(FA)  by  the  shaded  areas.  Now,  if  the  threshold  is  lowered,  as  in 
the  second  figure,  then  P(D)  shows  some  increase,  but  the  P(FA)  increases 
rather  drastically:  lowering  the  threshold  will  increase  the  probability  of 
a  detection,  but  the  "clutter"  will  also  increase.  On  the  other  hand,  if  the 
noise  curve  and  threshold  are  maintained  constant,  but  the  signal  gets  stronger, 
then  the  cxirve  of  signal  and  noise  will  move  to  the  right.  This  means  that  the 
detection  index  increases,  and  the  P(D)  increases  while  P(FA)  remains  the  same. 
This  says  no  more  than  that  a  stronger  signal  has  a  greater  probability  of 
being  detected. 

If  the  statistical  properties  of  the  noise  and  the  combined  signal  and  noise 
aire  known,  then  it  is  possible  to  make  curves  of  the  P(D)  and  P(FA)  relationships 
for  various  values  of  detection  index.  These  are  known  as  Receivers'  Operating 
Characteristics  (ROC)  curves. 
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P(D) 


RECEIVER  OPERATING  CHARACTERISITCS 
ROC 


P(FA) 

Typical  ROC  Curves 
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FALSE  ALARM  RATE 


The  rate  of  false  alarms  is  the  product  of  the  sampling  rate  and  P(FA). 
Since  sampling  rates  are  generally  high,  the  importance  of  a  very  low 
P(FA)  should  be  noted. 

For  example:  If  the  output  of  each  of  10  parallel  filters  is  sampled 
for  one  second  each,  a  P(FA)  of  1$  will  result  in  an  average  of  one 
false  alarm  every  10  seconds. 

The  sketch  below  shows  a. representative  plot  of  the  relationship  between 
d  and  P(FA),  for  a  P(D)  -  0.$; 


P(FA) 
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ENERGY  DETECTION 


For  a  passive  detection  system, 

a  completely  unknown  signal,  and  Gaussian  noise,  best  performance  is 
obtained  by  filtering  in  a  bandwidth  w,  taking  the  square  of  the  signal, 
and  accumulating  over  a  processing  time  t.  If  the  signal-to-noise 
ratio  is  small,  it  can  be  shown  that  the  detection  threshold  for  this 
process  is 


DTg  ■  5  lof  (d/wt) 


where  wt  is  called  the  time -bandwidth  product. 

To  decrease  DT  for  a  fixed  d,  it  is  necessary  to  increase  the  time -bandwidth 
product. 


Footnote  * 

The  detection  threshold  used  in  this  material  is  that  calculated  over 
the  complete  bandwidth  of  the  processor.  It  often  occurs  that  what  is 
specified  is  the  detection  threshold  in  a  one  Hertz  bandwidth,  Dr(l). 

This  latter  quantity  occasionally  leads  to  SONAR  equations  which  have  a 
slightly  sin5)ler  format  when  it  is  assumed  that  the  spectrum  levels  are 
constant  over  the  bandwidth  of  the  system.  Otherwise,  we  feel  that  the 
use  of  nr(l)  introduces  more  possibilities  for  confusion,  so  that  its  use 
creates  more  problems  than  it  solves.  (In  particular,  if  Urick's  excellent 
book  is  studied,  it  must  be  fully  understood  that  whereas  he  uses  DT  throughout 
most  of  the  text,  his  chapter  devoted  to  the  detection  threshold  uses  DT(lk).) 
The  formula  relating  the  two  quantities  is 
« 

DT(1)  -  DT  +  10  log  w 
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DOPPLER  SHIFT 


C3’ 


AJ  --  f 


While  decreasing  w  will  improve  the  ability  of  a  paasive  SONAR  to  detect 
tonals,  it  cannot  be  made  smaller  than  the  expected  Doppler  shifts. 

For  a  passive  SONAR  system,  this  requires  that  the  minimum  acceptable 
bemdwidth  be 

w(in  Hz)  =  0.7  X  (frequency  in  kHz)  x  (range  rate  in  knots) 
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EXAMPLE  1 


A  submarine  with  a  160  dB  re  l/<Pa  line  at  500  Hz  crosses  a  convergence 
zone  at  UO  km  from  an  omnidirectional  hydrophone.  The  bandwidth  of  the 
receiver  is  100  Hz  and  the  sea  state  is  3»  Assuming  a  convergence  gain 
of  12  dB,  find  out  how  long  the  signal  must  be  accumulated  to  give  a 
P(D)  -  0.50  and  a  P(FA)  -  10“^. 


Solution: 

TL  «  20  log  r  +  ar  -  0  =  20  log  (U  x  IcA)  ♦  5  x  10  ‘LX  ic/^  -  12  *  82  d£ 

ML  •  NSL  +  10  log  w  -  DI  -  66+20-0  -  86  dB 
SL  -  TL  >  ML  +  DT 

DT  SL  -  TL  -  ML  =  160  -  82  -  86  -  -8  dB 

DT  ■  5  log  (d/wt) 
where  d  •  l5 

log  (wt/d)  ■  1,6 
wt/d  ■  40 
wt  •  600 
t  -  6  sec 
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EXAMPLE  2 


Sams  as  Example  but  let  w  be  just  vide  enou|,h  to  encompass  a  20  knt 
target. 

Solution: 


where 


w 

- 

0.7  X  (frequency  in  kHz)  x  (range  rate  in 

knots) 

- 

0.7  k  20  -  7  Hz 

ML 

3 

NSL  +  10  log  w  -  DI  -  66+8.5- 

74.5  dB 

DT 

< 

SL  -  TL  -  ML  -  160  -  82  -  74.5  - 

3.5  dB 

DT 

- 

5  log  (d/wt) 

d 

m 

15 

log  (d/wt)  »  0.7 


d/wt  ■  5 

wt  -  3 
t  ■  3/7  sec 
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EXAMPLE  3 


A  conventional  submarine  at  periscope  depth  is  traveling  at  4  kts.  The 
Sea  State  is  3  and  the  layer  depth  is  100  m.  At  what  range  will  the 
submarine  be  detected  by  a  36  m  deep  hydrophone  listening  at  1000  Hz 
if  DI  =  20  dB  and  DT  =  0  dB? 

Solution: 

SSL  -  TL  ^  NSL  -  DI  +  DT 

TL  ^  SSL  -  NSL  +  DI  -  DT  =  120  -  62  +  20  =  78  dB 

For  sound  to  be  trapped  in  the  duct, 

f  >  f  .  =  2*10^A>^'^^  =  200  Hz 

^  min 

So  we  use  the  TL  for  a  mixed  layer 


TL  = 

10  log  + 

10  log  r  +  ar  +  b(r/r  ) 

s 

ri¬ 

ll 

1.05-10^  D/(d 

1/2  2 

-  z  )  '  =  1.05-10  (100/8)  =  1310  m 

s 

r  = 

8.4-10^ 

=  8. 4 *10^ *10  =  8400  m 

s 

From  the  graphs  and  equations 
a  =  6*10  ^  dB/m 

b  =  1.7  dB /bounce 

Notice  that  ar  +  b(r/r  )  +  (a  +  b/r  )  r 

s  s 

Thus  78  =  10  log  1310  +  10  log  r  + 

10  log  r  +  2.6*lo"^  r  =  47 

Now,  trial  and  error 


r 

10  log  r 

2.6‘10”^  r 

TL 

5*10'* 

47 

13 

60 

2-10^ 

43 

5.2 

48 

4 

1.5-10 

41.8 

3.9 

46 

4 

1.7-10 

42.3 

4.4 

1 

47 

r  =  1.7-10'^  m 


=  2.6*10”^  r 

2.6-10"'^  r 
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EXAMPLB  h 


Repeat  Example  3  with  DI  ■  0  dB. 

Solution: 

TL  ^  SSL  -  NSL  +  DI  -  DT  -  120  -  62  •  58  dB 

-4 

TL  •  10  log  r  +  10  log  r  +  2.6  x  10  r 
58  -  10  log  1310  +  10  log  r  +  2.6  x  lo"^  r 


58  •  31  •» 

10  log  r  + 

2.6  X  lo"^ 

r 

10  log  r  + 

2.6  X  lO"'^  r 

=  27 

r 

10  log  r 

2.6  X  lO"'^ 

TL 

5  X  10^ 

27 

0.13 

27 

But  r  < 

So  we  must  use  this  form  of  the  TL.' 

TL  »  20  log  r  ♦  ar 

58  -  20  log  r 

log  r  -  2.9 

r  •  800  m 


Bingo 
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EXAMPLE  5 


Repeat  Example  3  with  f  ■  100  Hz.  Assume  "heavy  shipping" as  would  be 
encountered  in  the  Mediterranean  Sea. 


Solution: 


f  <  f , 


min 


No  trapping  in  the  layer.  Use 

TL  »  20  log  r  +  ar 

TL  SSL  -  NSL  +  DI  -  DT  •  133  -  76  +  20  -  77  dB 
20  log  r  +  10**^  r  -  77 


r 

20  log  r 

10"^  r 

TL 

8  X  10^ 

78 

8  X  10"^ 

78 

7  X  10^ 

77 

7  X  10'^ 

77 

r  •  7  X  10^  m 


Bingo 
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ACTIVE  SONAR 


SL  -  2TL  +  TS  ML  +  DT 

n 


vhere 

SL  ■  source  level  of  transmitter 
21'L  •  tvo-vay  transmission  loss 
TS  ■  target  strength 


Sacb  term  in  the  above  equation  vill  now  be  discussed  separately. 


Notet  the  above  active  SONAR  equation  is  for  a  monostatic  case, 
wherein  the  source  and  receiver  are  coincident  so  that  the  TL 
from  source  to  target  is  the  same  as  the  TL  from  target  to 
receiver.  Our  discussion  will  usually  be  restricted  to  this 
case.  Extension  to  the  bistatic  situation  in  which  source  and 
receiver  are  at  different  locations  is  trivial:  If  the 
transmission  loss  from  source  to  target  is  TL,  and  from  target 
to  receiver  TL' ,  then  all  that  is  required  is  to  replace 
2TL  with  TL  +  TL’ ,  so  that  we  have 


SL  -  TL  -  TL’  +  TS  ^  ML  +  DT 
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LIST  OF  ADDITIONAL  SYMBOLS  FOR  ACTIVE  SONAR 


SL  *=  source  level  of  the  active  sonar  system 
TS  ■  target  strength  of  the  desired  target 
RL  =  reverberation  level 

TS'  *  target  strength  of  the  scatterers  causing  reverberation 
TL  =  transmission  loss  from  source  to  target 

TL'  =  transmission  loss  from  target  to  receiver,  equal  to  TL  in 
the  case  of  monostatic  systems 


Lg  =  "effective  length"  of  a  target 

t  *  pulse  duration  of  the  acoustical  signal  from  the  source 
s^  =  a  measure  of  the  scattering  cross-section  per  unit  volume 
Sg  =  a  measure  of  the  scattering  cross-section  per  unit  area 
Sy  *  scattering  strength  for  volume  scatterers 
S_  »  scattering  strength  for  surface  scatterers- 
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(§)  -  2TL  +  TS  ^  ML  +  DT 
SOURCE  LEVEL  -  ACTIVE  SONAR 


SL  ■  sound  pressure  level  1  m  from  the  acoustic  center  of  the 
transmitting  transducer. 


SL  ■  10  log  (Pover)  +  Iff  +  171 


where 


or 


where 


and 


Power  •  Acoustic  power  output  of  the  transducer  in  Watts 

DI  =  10  log  (the  ratio  of  the  intensity  on  the 

acoTistic  axis  to  the  intensity  at  the  same 
distance  from  an  omnidirectional  transmitter 
with  the  same  acoustic  power  output) . 


m 

^o(r) 


20  log  P(r)/Po(r) 

pressure  produced  by  an  (xnnidirectlonal  transducer 
with  the  same  power  output 


171  ■  factor  to  take  care  of  the  units 


cd 

Oh 

210 

K 

r-l 

S 

220 

% 

200 

O 

Q 

180 

1 

CO 

160 

Acoustic  power  output  (Watts) 
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SL  -  2TL  +  ML  +  DT 

TARGET  STRENGTH 


TS  •  (sound  pressure  level  of  reflected  wave  1  m  from  the 
acoustic  cen^^r  of  the  target  in  a  direction  back 
towards  the  receiver)  -  (sound  pressure  level  of 
the  incident  wave  at  the  position  of  the  target  if 
the  target  were  absent). 

or 

TS  •  20  log  P*(r'"l)/P!^r) 

where 

P'(r’«l)  is  the  pressure  of  the  reflected  wave  1  m  from  the 
acoustic  center  of  the  target  in  a  direction  back  towards  the 
receiver. 


Target  Strength  depends  on  shape,  construction,  and  orientation  of  the 
target,  and  on  the  frequency  and  pulse  length  of  the  sound  wave. 
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EFFECT  OF  TARGET  SIZE 


Because  of  the  effective  length  L  of  the  target,  portions  of  the 
incident  pulse  which  reflect  from^the  tail  (in  the  above  sketches) 
travel  a  total  distance  of  2L  further  than  the  reflections  from 
the  nose.  This  means  that  thl  apparent  duration  of  the  pulse 
is^  stretched  from  t  to  (t  +  2L  /c)  .  For  very  long  pulses, 
the  additional  time  interval  cofresponding  to  the  round  trip, 

2L  /c,  can  be  neglected.  For  short  pulses,  however,  this  stretching 
can  be  quite  important  in  determining  the  target  strength. 
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The  Table  shovs  a  few  representative,  but  very  approximate,  values  of 
target  strength  for  a  variety  of  objects. 


Target 

Frequency 

in  kHz 

Fleet  Submarine 

10 

(approx  100  m  long) 

2k 

S  Class  Submarine 
(approx  70  m  long) 

2k 

ToiTpedo 

2k 

Surface  ''essel 

3  ft  mine 

30-90 

Approximate  Target  Strength  (dB) 
Beam  Aspect  Bow-Stern  Aspect 

35  30 

2$  16 

20  13 

-20 

10  to  30 
-8 


Extracted  from:  Physics  of  Sound  in  the  Sea,  NDRC 

Principles  of  Underwater  Sound  for  Engineers,  Urick 
Principle  and  Applications  of  Underwater  Sound,  NIEC 
NAVSHIPS  0967-129-3010 
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SL 


2TL  +  TS  ^  +  DT 


MASKING  LEVEL 

The  undesired  noise  competing  with  the  desired  signal  can 
arise  from  two  possible  mechanisms  in  the  case  of  active  SONAR t 

Masking  from  ambient  or  self  noise 

Masking  from  reverberation 

Which  of  these  two  mechanisms  is  the  dominant  depends  on  the 
range  between  source  and  target.  Generally  speaking,  for  low- 
powered  active  systems  the  performance  will  be  noise-limited, 
in  that  the  maximum  effective  detection  range  will  be  determined 
when  the  echo  level  diminishes  below  that  level  for  which  it  can 
be  extracted  from  the  ambient  or  self  noise.  For  higher-powered 
active  systems,  however,  reverberation  becomes  an  important  source 
of  interference.  Reverberation  tends  to  decrease  with  increasing 
range,  but  less  rapidly  than  the  echo  level.  Thus,  if  the  echo 
level  has  diminished  until  it  is  buried  in  the  reverberation, 
but  the  reverberation  has  not  yet  fallen  sufficiently  below 
the  ambient  (or  self)  noise,  then  the  maximum  detection  rangp 
is  determined  by  the  reverberation.  The  two  situations  are 
suggested  in  the  figures: 


The  choice  is  a  little  more  subtle  than  suggested  by  the  simple 
curves  above.  The  detection  threshold  for  noise  masking  may 
be  different  than  that  for  reverberation  masking,  Ahd  the 
detection  threshold  for  reverberation-limited  performance  may 
depend  strongly  on  the  doppler-induced  difference  between  the 
frequency  of  the  reverberant  sound  and  that  reflected  from  the 
desired  target. 
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AMBIENT  AND  SELF  NOISE 


The  masking  level  is  given  by 

ML  =  NL  -  AG 

where  the  noise  level  NL  is  that  from  the  ambient  noise  of  the 
ocean,  noise  generated  by  the  receiving  platform,  or  the  combination 
of  the  two  if  they  are  of  nearly  equal  strengths. 

Recall  that  if  the  noise  spectrum  level  is  sufficiently 
smooth  over  the  receiver  bandwidth,  then  we  can  write 

NL  =  NSL  +  10  log  w  . 

Otherwise,  the  total  bandwidth  must  be  broken  down  into  smaller 
bands,  the  band  level  found  for  each,  and  then  these  band  levels 
combined  according  to  the  nomogram  for  combining  independent 
levels. 
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REVERBERATION 


If  the  sound  field  generated  by  the  source  is  strong  enough,  and 
there  are  enough  scatterers  present  in  the  vicinity  of  the  desired 
target,  then  unwanted  reflections  and  scattering  of  sound  from 
the  undesired  scatterers  are  received  along  with  the  echo  from 
the  desired  target.  The  echo  level  of  the  signal  received  from 
the  desired  target  is 

EL  *  SL  -  2TL  +  TS  (monostatic) 

where  TS  is  the  target  strength  of  the  desired  target.  The 
reverberation  level  RL  can  be  written  in  the  form 

RL  *  SL  -  2TL  +  TS’  (monostatic) 


where  TS’  is  the  combined  target  strength  of  all  the  undesired 
targets  which  are  generating  the  masking  level.  Since  the  RL 
is  the  masking  level  when  reverberation  is  the  dominant  source 
of  unwanted  sound,  the  SONAR  equation  becomes 


or 


EL  -  RL  ^  DT 
TS  -  TS’  DT 


(monostatic ) 


The  central  problem  in  obtaining  TS’  is  determining  the  volume 
or  surface  area  containing  the  scatterers  which  is  illuminated 
by  the  source,  sending  sound  to  the  receiver,  and  sending  the 
sound  so  that  it  arrives  at  the  same  time  as  the  signal  from  the 
desired  target.  This  is  suggested  in  the  sketch  below. 
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VOLUME  REVERBERATION 


The  target  strength  TS’  of  the  undesired  scatterers  can  be 
written  as 

TS’  *  10  log  s^V 

where  V  is  the  volume  within  which  the  reverberation  sent  to 
the  receiver  is  in  competition  with  the  reflected  signal  from 
the  desired  target  and  Sy  is  a  measure  of  the  scattering  cross 
section  per  unit  volume.  It  has  become  conventional  to  express 
Sy  in  dB  terms  through  the  definition 

Sv  =  10  log  Sv 

where  Sy  is  the  scattering  strength.  Thus,  we  have 

TS’  *  Sy  +  10  log  V. 

The  calculation  of  the  reverberating  volume  V  contains  the  heart 
of  the  problem.  As  suggested  previously,  this  is  a  study  in 
geometry,  involving  the  way  in  which  major  lobes  of  source  and 
receiver  overlap. 

If  we  restrict  ourselves  to  the  monostatic  situation,  and  assume 
that  the  major  lobe  of  the  receiver  completely  overlaps  the 
major  lobe  of  the  source,  it  then  follows  that  the  volume  V  is 
given  by 

10  log  V  *  10  log  Ct  +  +  ^0  log  SI 

where  t  is  the  time  duration  of  the  acoustical  pulse  and  XL  is 
the  solid  angle  subtended  by  the  major  lobe  of  the  radiation 
pattern  of  the  source.  Recall  that  this  is  related  to  the 
directivity  index  by  DI  10  log(4Tr/fV) .  The  term  TLg  arises 
from  the  fact  that  the  cross  sectional  area  of  the  reverberating 
volume  clearly  depends  on  the  range  from  the  source,  and  this  can 
be  seen  to  be  expressed  by  the  geometrical  spreading  of  the  sound 
beam.  Thus,  for  this  simple  monostatic  situation,  we  have 

TS’  «  Sy  +  10  log  ct  +  TLg  +  10  log  kw  -  DI 
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SURFACE  REVERBERATION 


In  the  case  of  scattering  which  arises  from  a  surface  (such  as 
the  sea  surface  or  the  bottom)  or  from  some  well-defined  layer 
such  as  the  bubble  layer  underlying  the  water-air  interface, 
the  reverberation  target  strength  is  expressed  as 

TS’  *  Sg  +  10  log  A 

where  A  is  the  surface  area  generating  reverberation  in  competition 
with  the  desired  target  and  Sg  is  the  scattering  strength  (per 
unit  area). 

As  before,  the  problem  is  one  primarily  of  geometry,  involving 
the  way  in  which  source  and  receiver  major  lobes  overlap.  If  we 
again  restrict  ourselves  to  the  simple  monostatic  geometry  used 
previously  (for  which  the  major  lobe  of  the  receiver  completely 
overlaps  that  of  the  source),  then  the  surface  area  of  interest 
can  be  seen  to  be 


A  =  i  r  0  c  t 

where  0  is  the  angle  in  radians  of  the  horizontal  angular  width 
of  the  radiation  pattern  of  the  source.  This  then  gives  the 
equation 

TS*  =  +  10  log  ^  +  10  log  r  +  10  log  0 

s  2 
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REPRESENTATIVE  SCATTERING  STRENGTHS 


The  values  quoted  here  are  representative  only.  Large 
deviations  can  occur  depending  on  location,  season,  biological 
activity,  and  so  forth. 

Deep  Scattering  Layers:  The  scattering  strength  S^ 
can  vary  between  about  -90  to  -60  dB  within  the 
various  sublayers. 

Water  Volume:  In  the  clear  open  ocean  waters,  Sy 
can  be  found  to  range  between  -100  to  -?0  dB. 

Surface  Scattering:  For  reasonably  low  grazing  angles, 
the  combined  scattering  from  the  ocean  surface  and 
the  population  of  bubbles  found  at  shallow  depths 
below  the  surface  gives  scattering  strengths  Sg 
ranging  around  -50  to  -30  dB  for  sea  states 
between  1  and  4. 

Bottom  Scattering:  For  grazing  angles  between  20° 

and  60°,  representative  values  of  the  scattering 
strength  Sg  range  between  -40  and  -10  dB.  The 
exact  values  depend  very  strongly  on  the  particular 
bottom  and  its  composition.  In  general,  however, 
as  the  grazing  angle  goes  to  zero,  the  scattering 
becomes  negligible  and  Sg  goes  to  - co  dB. 


REVERBERATION  LEVEL  DEPENDENCIES 

Notice  that  for  either  volume  or  surface  reverberation 
the  reverberation  target  strength  TS'  increases  with  increasing 
range,  either  as  TL^  (volume)  or  10  log  r  (surface).  This  is  to 
be  compared  with  ambient  noise  masking,  for  which  there  is  no 
range  dependence.  This  means  that,  given  sufficient  source  level, 
as  increases  SONAR  performance  will  change  from  ambient- 

noise  limited  performance  to  reverberation  limited  performance. 
While  detection  for  ambient  noise  limited  performance  depends 
on  the  source  level  and  directivity  index,  notice  that  when  the 
performance  becomes  reverberation  limited  there  is  no  dependence 
on  the  source  level:  additional  power  will  not  improve 
reverberation  limited  performance.  Reverberation  can  be  reduced 
by  increasing  the  directivity  of  the  source  and/or  receiver,  or 
by  decreasing  the  pulse  length,  given  by  ct.  Increased  directivity 
has  the  disadvantage  of  reducing  scanning  rates,  and  decreased 
pulse  length  may  result  in  reduction  of  the  target  strength  TS 
of  the  desired  target. 
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DOPPLER  SHIFT  FOR  ACTIVE  SONAR 


In  the  case  of  active  SONAR,  any  motion  of  the  platform 
containing  the  source  will  cause  the  signal  propagated  into 
the  water  to  have  a  frequency  shift  which  depends  on  the 
bearing  into  which  the  signal  is  sent.  This  doppler-shifted 
signal  is  subjected  to  a  further  doppler-shift  when  it  is 
reflected  from  the  moving  target.  Without  going  into  the  specifics 
of  the  various  frequency  shifts,  we  shall  simply  quote  some  results 
of  interest: 

Because  of  the  two-way  propagation  of  acoustical  energy  from 
the  moving  source  to  the  moving  target  and  then  reflection  back 
to  the  moving  receiver,  the  required  bandwidth  for  the  receiver 
of  an  active  sonar  operated  under  monostatic  conditions  is 

w  in  Hz  =  1.4* (Range  rate  in  knots )• (Frequency  in  kHz) 

Notice  that  this  is  twice  the  bandwidth  required  in  the  passive 
SONAR  case. 

If  the  target  is  stationary  with  respect  to  the  surrounding 
water,  then  the  desired  echo  will  have  the  same  frequency  as  the 
reverberation  and  the  ability  of  the  receiving  system  to  recover 
signal  from  noise  will  be  minimized.  On  the  other  hand,  when  the 
target  is  moving  through  the  water,  then  the  echo  and  reverberation 
will  have  different  frequencies  and  detection  will  be  much  easier. 
This  effect  is  quite  easy  to  observe  by  eari  The  detectibility 
of  an  "up-doppler"  or  "down-doppler"  echo  in  reverberation  is 
very  easy  even  when  the  relative  loudness  of  the  desired  echo  is 
rather  weak.  When  there  is  no  frequency  change  between  echo  and 
reverberation,  however,  the  ability  of  the  ear  to  detect  the 
echo  is  drastically  reduced.  The  strong  increase  of  the  detection 
threshold  for  slow  targets  points  to  the  use  of  parallel  filters 
of  relatively  narrow  bandwidths  to  facilitate  the  separation  of 
signal  from  reverberation  when  the  difference  in  frequency  is 
small. 
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SL 


2TL  +  TS  ^  ML  + 


DETECTION  THRESHOLD 


Energy  Detection 

If,  as  was  the  case  for  the  passive  SONAR  equations,  the 
mechanism  for  the  detection  of  the  desired  signal  is  based  on 
detecting  an  excess  of  acoustical  energy  in  the  frequency  band 
of  interest,  then  we  have  the  same  expression  for  the  detection 
threshold , 


DT  =  5  log  ^ 

where  t  is  the  integration  time.  In  the  case  of  active  sonar, 
t  should  be  the  length  of  the  received  echo.  If  the  integration 
time  is  less  than  or  greater  than  this  pulse  duration,  it  can  be 
shown  that  the  DT  is  increased  above  this  expression,  and 
performance  degraded. 

Correlation  Detection 


An  alternative  mode  of  signal  processing  is  possible  in 
the  case  of  active  SONAR.  Since  the  amplitude  and  frequency 
properties  of  the  tone  burst  generated  by  the  source  are  known, 
it  is  possible  to  search  for  a  signal  of  these  same  properties 
in  the  received  echo.  Ifl  "the  detailed  shape  of  the  received 
echo  matches  that  of  the  sent  pulse,  then  it  can  be  shown  that 
the  detection  threshold  is  given  by 

DT  =  10  log  ^ 

when  the  processing  is  that  of  correlation  of  the  received  signal 
against  a  replica  of  the  generated  pulse.  There  are  many 
modifications  of  this  basic  idea,  but  all  rely  on  the  technique 
of  multiplying  the  received  signal  by  a  time-delayed  model  of 
the  sent  pulse  and  integrating  the  resultant  product  over  the 
pulse  duration.  This  is  done  for  monotonically  increasing 
delay  times,  and  the  resulting  function  of  delay  time  studied 
for  the  presence  of  a  sharp  peak  representing  the  matching  of 
the  delayed  replica  with  the  same  signal  appearing  in  the  received 
signal  and  noise.  In  practice,  the  finite  size  of  the  target, 
multipath  interference,  and  the  fluctuations  introduced  by  the 
inhomogeneities  of  the  transmission  of  the  signal  through  the 
water  alter  the  received  echo  so  that  the  correlation  process 
is  not  as  good  as  would  be  desired.  As  a  result,  the  DT  for 
a  real  system  operating  in  the  real  world  is  more  than  the 
theoretical  expression  would  predict. 
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SUMMARY  OF  THE  ACTIVE  SONAR  EQUATIONS 


Ambient  noise  limited  conditions 


SL  -  2TL  +  TS  ^  NBL  -  AG  +  10  log  w  +  DT 


Reverberation  Umited  conditions 


Ts  -  TS*^  nr 


107 


EXAMPLE 


An  active  SONAR  operates  at  1  kHz  with  a  source  level  of  220  dB  re  1  >^Pa, 
a  directivity  index  of  20  dB,  a  horizontal  beam  width  of  10°,  and  a  pulse 
length  of  0.1  sec.  Correlation  detection  is  used  and  it  is  desired  to  have 
P(D)  -  0,50  with  P(FA)  »  10~ii,  The  target  strength  of  the  submarine  is 
expected  to  be  30  dB  and  its  speed  may  be  up  to  20  kts.  Both  the  SONAR 
and  the  submarine  are  in  a  mixed  layer  of  depth  100  m  with  the  source  at 
36  m.  The  sea  state  is  3>  the  scattering  strength  for  surface  scatterers 
is  -30  dB,  and  volimia  scattering  is  negligible.  For  simplicity,  assume  that 
conditions  are  such  that  the  detection  threshold  is  the  scime  for  both  noise- 
limited  and  reverberation-limited  performance. 

Solution: 

w  -  l.ii  (frequency  in  kHz)  (range  rate  in  kts)  -  I.I4  x  1  x  20  •  28  Hz 

d  ■  15  same  calculation  as  in  example  1  (passive) 

DT  •  10  log  (d/2wt)  -  10  log  (15/2x28x0.1)  -  10  log  2.68  -  5  dB 

NL  -  NSL  -  DI  +  10  log  w  -  62  -20  +  lit  -  56  dB 

r^  ■  1310  m  and  r^  ■  84OO  m  same  calculation  as  Example  3  (passive 

Assume  that  the  detection  is  noise  limited 

2TL  SL  +  TS  -  NL  -  nr  -  220  +  30  -56  -  5  -  I89 

TL  95 

-4 

TL  •  31  +  10  log  r  +2.6  X  10  r  same  calculation  as  Example  3 

10  log  r  +  2.6  X  10  ^  r  =  95-31  =  64 


r 

10  log  r 

2.6  X  10”^  r 

TL 

10^ 

40 

26 

43 

10^ 

50 

26 

76 

5  X  lo"^ 

47 

13 

60 

6  X  10^ 

48 

16 

64 

r  =  6  X  10^  m  if  noise  limited 
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EXAMPLE  (CONTINUED) 


But,  what  is  the  RL  at  this  range? 

TS'  •  Sg  +  10  log  (ct/2)  +  10  log  r  +  10  log  G 

»  -30  ♦  10  log  (1500  X  0.1/2)  +  48  +  10  log  O.I7I4 

■  -30  +19+48-8  ■  29  dB 

RL  -  SL  -  2TL  +  TS ' 

■  220  -  189  +  29  .  60  dB 
So 

RL  NL 

Therefore  the  detection  will  be  reverberation  limited 
TS  -  TS  •  ^  DT 

TS  •  •  -30  +  19  +  10  log  r  -  8  -  -19  +  10  log  r 

TS  -  TS'  >  ET 

30+19  -  10  log  r  5 

10  log  r  ^  iiii. 

r  -  2.5  X  10^  m 
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APPENDIX  A:  LOGS 


All  you  have  ever  wanted  to  know  about  the  decibel  but  were  afraid  to  ask. 
The  language  in  which  the  SONAR  equation  is  written  is  the  language  of  the 
/'  decibel.  This  self-paced  exercise  is  intended  to  make  you  familiar  with  the 

grammar  of  this  language. 

Instructions .  For  each  numbered  section,  first  attempt  to  do  the  problems. 

(The  answers  are  on  the  last  page.)  If  you  have  no  troubles,  skip  over  to  the 
next  section.  On  the  other  hand,  if  you  have  difficulties^  read  through  the 
discussion  and  then  try  the  problems  again.  Go  to  the  next  section  only  when 
you  can  do  the  problems  without  difficulty.  By  the  time  you  finish  all  the 
sections,  you  should  be  speaking  "decibel"  like  a  native. 

1.  EXPONENTIAL  NOTATION 

Problems  A 

1.  Write  each  of  the  following  numbers  in  exponential  forms 

a)  2,365,000 

b)  872 

c)  3340 

d)  42.1 

e)  600 

f)  1.0 

2.  Write  each  of  the  following  numbers  in  conventional  form 

a)  2.689  X  10^ 

b)  7.30  X  10^ 

c)  6.29  X  lo"* 

d)  8.9  X  10^ 

e)  3.261  X  10^ 

f)  1.00  X  10° 


no 


3.  Write  each  of  the  following  numbers  in  exponential  form 

a)  0.00252 

b)  0.012 

c)  0.00002 

d)  0.01002 

e)  0.00328 

4.  Write  each  of  the  following  nianbers  in  conventional  form 

a)  1.03  X  lo"^ 

b)  3.2  X  10~^ 

c)  6.89  X  lo"'* 

d)  1.003  X  lo"^ 

e)  1.00  X  10“^ 


A  succinct  way  of  expressing  very  large  and  very  small  numbers  is  by 
exponential  (or  scientific)  notation.  For  exaimpler  one-hundred  and  twenty-six 
million  can  be  written  either  as 

126,000,000 

or  as 

1.26  X  10®  , 

where  the  exponent  8  denotes  the  number  of  places  that  the  decimal  must  be 
moved  to  the  right  to  express  the  nvimber  in  its  conventional  form. 

There  is  no  unique  way  of  expressing  a  nxmber  in  exponential  form,  since 
1.26  X  10®  =  12.6  X  lo’  =  0.126  x  10®  =  etc. 

Which  form  you  use  depends  on  convenience. 

Numbers  smaller  than  unity  can  be  handled  in  a  similar  manner.  For 
example,  one-hundred  and  twenty  millionths  can  be  written  either  as 


0.000120 


where  the  negative  exponent  denotes  that  the  decimal  place  is  to  be  moved  to 


the  left. 

A  special  case  that  must  be  watched  is 

1  =  1  X  10°  =  10° 

which  is  consistent  with  the  above  conventions. 
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2 .  ALGEBRA  OF  EXPONEMTS 


Problems  B 

1.  Calculate  the  following  and  express  the  answer  in  exponential  form 

a)  6.15  X  10^  +  2.34  X  10^ 

b)  3.2  X  10^  -  1.46  X  10^ 

c)  9.1  X  10^  +  2.4  X  10^ 

d)  3.58  X  lo”^  +  1.26  x  lo"^ 

e)  6  X  lo"^  +  3.281  x  lo”'^ 

f)  3.5  X  lo'*  -  2  X  10^ 

2.  Calculate  the  following  and  express  the  answers  in  exponential  form 


a) 

2.50 

X 

10^  X  2.00  X  lo'* 

b) 

4.40 

X 

10^  /  2.20  X  10^ 

c) 

3.26 

X 

10^  X  2.00  X  lo”^ 

d) 

5.34 

X 

lo"^  /  2.67  X  10^ 

e) 

8.68 

X 

lO”^  /  4.34  X  lo”^ 

3.  Calculate  the  following  and  express  the  answer  in  exponential  form 

a)  (1.1  X  lo'*)^ 

b)  (2.0  X  lo”^)'* 

c)  (1.44  X  10^) 

-2  1/2 

d)  (1.44  X  10 

e)  (8.0  X  10^)^^^ 

f)  (2.56  X  10^)°*^ 

t  /«  0.3333 _ 

g)  (8  X  10  ) 

To  add  and  subtract  nvimbers  expressed  in  exponential  form  it  is  necessary 

4 

to  express  both  numbers  in  the  same  exponent.  For  example,  adding  3.56  x  10 

3 

to  2.1  X  10  IS  accomplished  by 

3.56  X  10^ 

0.21  X  10^ 

3.77  X  lo"* 
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When  asked  to  subtract  8.67  x  10  ^  from  2.23  x  10^,  note  that  the  answer  is 
2.23  X  10^,  for  the  term  to  be  subtracted  affects  the  larger  term  beyond  the 
figures  retained. 

Exponential  notation  really  shows  its  worth  in  calculations  involving 
multiplication  and  division.  We  see  that 

120  X  2000  =  240,000 

or  1.2  X  10^  X  2  X  10^  =  2.4  x  10^ 

and 

36000  /  200  =  180 
or 

3.64  X  10^  /  2  X  10^  =  1.8  x  10^. 

The  rules  are  simple: 

In  multiplication,  exponents  add. 

In  division,  exponents  subtract. 

A  special  case  that  is  often  useful  is 

2  0  2  -2 
1  /  10  =  10  /  10^  =  10  , 

which  allows  exponents  to  be  moved  between  numerator  and  denominator. 

Power  and  roots  are  trivial  in  exponential  notation: 

(1.2  X  10^)^  =  1.2  X  10^  X  1.2  X  10^  =  1.44  x  10®, 


or  for  any  power  n 


and. 


or  for  any  power  of  n 


do”')^  =  lo"^ 


(1.44  X  =  (1.44)*''=  X  (10")'''^  ■  ' 


This  latter  form  introduces  the  possibility  of  fractional  exponents;  these  will 
be  discussed  in  Section  4. 


I 

I 

1 
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3.  LOGS  (IN  BASE  TEN) 

Problems  C 


Find  X 

a) 

log 

10^  =  X 
-6 

b) 

log 

10  =  X 

c) 

log 

10^  =  4 

d) 

log 

10*  =  -3 

e) 

log 

10*  =  0 

Since  we  say  that  "ten  to  the  fourth  power  equals  ten  thousand"  we  can 
ask  "to  what  power  must  ten  be  raised  to  get  ten  thousand?"  The  answer  is 
obviously  "four."  You  are  now  an  expert  on  logarithms,  since  the  "log"  of 
10,000  is  4.  This  is  written 

log  10,000  =  4  or  log  lo"*  =  4 
and  this  is  all  there  is  to  logs. 

Find. the  log  of  x  by  asking  yourself  to  what  power  must  I  raise  10  to 
get  X. 

It  is  equally  appropriate  to  ask  "the  log  of  what  is  4,"  that  is  what 
is  X  so  that  log  10  =4. 
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4.  DECIMAL  EXPONENTS 


PROBLEMS  D 

1.  Find 

a)  log  4 

b)  log  40 

c)  log  0.4 

d)  log  4^ 

e)  log  1/4 

2.  Find  the  number  whose  log  is 

a)  0.70 

b)  1.70 

c)  -1.30 

d)  -1.70 

e)  -2.70 

So  far  we  have  assumed  integer  exponents,  but  everything  we  have 
said  also  is  true  for  decimal  exponents.  For  example, 

10"-^"=  2.0 

and 


log  2  =  0.30 

The  table  below  will  allow  you  to  calculate,  with  a  fair  degree  of  accuracy, 
all  logs. 


X 

1.0 

1.2 

1.5 

1.7 

2.0 

.3.0 

4.0 

5.0 

6.0 

7.0 

8.0 

9.0 

10.0 

lo^  X 

0.00 

0.08 

0.18 

0.23 

0.30 

0.48 

0.60 

0.70 

0.78 

0.85 

0.90 

0.95 

1.0( 
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If  you  ever  need  more  accuracy  you  will  have  to  use  slide  rule,  calculator, 
or  tables, 

To  find  the  log  of  a  number  greater  than  10  proceed  as  follows: 
log  7000  =  log  7  X  10^  =  log  7  +  log  10^ 

=  0.85  +3  =  3.85 

or 

log  0,003=  log  3  x  10  ^=  log  3  +  log  10  ^ 

=  0,48  -  3  =  -2.52 

Uiere  is  one  property  of  logs  that,  if  used  correctly,  will  greatly 

simplify  many  calculations.  This  will  be  illustrated  by  an  example: 

2 

log  36  =  log  6  =  log  (6x6)  =  log  6  +  log  6  =  2  log  6=  1.56 

f Calculate  log  36  and  compare  to  2  log  6.) 

In  general 

log  X  =  n  xog  x. 

It  is  also  true  that 

log  xy  =  log  x  +  log  y, 

log  36  =  log  (4x9)  =  log  4+2  log  3  =  0,60  +  0,96  =  1.56 
But,  remember  that 

log  (x  +  y)  is  not  equal  to  log  x  +  log  y  III 


117 


5.  THE  DECIBEL 


PROBLEMS  E 


1.  E:q>ress  these  signal- to-noise  ratios  in  dB 

a)  2 

b)  1/2 

c)  10 

d)  1/10 

e)  100 

2.  The  voltage  of  a  signal  is  1  volt*,  express  the  signal  level  with 
reference  to  the  following  voltages 

a)  lO"^^ 

b)  10"® 

c)  1 

d)  10® 


By  convention,  when  the  log  of  a  pressure  or  voltage  is  multiplied 
by  20,  the  resulting  number  is  said  to  be  measured  in  "decibels”,  donated 
by  dB. 

20  log  2  =  6dB 
20  log  10  =  20dB 

Decibels  are  used  in  two  different  ways: 

a)  to  express  a  ratio  such  as  the  signal-to-noise  ratio 

signal- to-noise  ratio  =  20  log  (signal  voltage/noise  voltage) 

b)  to  express  a  level  such  as  signal  level 

signal  level  =  20  log  (signal  voltage/reference  voltage) 

In  the  latter  case  it  is  absolutely  essential  to  specify  the  reference, 
e.g.,  "the  signal  level  is  120dB  re  10  ^  volts". 
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Problems  A: 


Problems  B: 


Problems  C: 


Problems  D: 


Problems  E: 


ANSWERS 

1. a)  2.365  X  10^,  b)  8.72  x  10^,  c)  3.34  x  10^,  d)  4.21  x  10^, 

e)  6  X  10^,  f)  1.0  X  10° 

2.  a)  268.9,  b)  730,  c)  62,900,  d)  8,900,000,  e)  326.1 

f)  1.00 

3.  a)  2.52  X  lo"^,  b)  1.2  x  lO"^,  c) 2  x  lo"^,  d)  1.002  x  lo"^ 
e)  3.28  X  10"°. 

4.  a)  0.00000103,  b)  0.032,  c)  0.000689,  d)  0.1003,  e)  0.100 


l.a) 

8.49  X  10^,  b)  -1.14  X  10^, 

c)  3.3  X  10°,  d)  4.84 

X  lO”^ 

e) 

3.34  X  lO"'^,  f)  3.5  X  lo"^. 

2. a) 

5.00  X  10^°,  b)  2.00  lO"^, 

c)  6.52  X  10^,  d)  2.00 

X  lO"^, 

e) 

2.00  X  10. 

3. a) 

1.2  X  10°,  b)  1.6  X  lO”^^, 

c)  1.2  X  10^,  d)  1.2  X 

lo"^. 

e) 

2.0  X  10^,  f)  1.6  X  10^,  g) 

2.  X  10‘^. 

l.a)  3,  b)  -6  ,c)  4,  d)  -3,  e)  0 

1. a)  0.60,  b)  1.60,  c)  -0.40,  d)  1.20,  e)  -0.60. 

2.  a)  5.0,  b)  50,  c)  0.050,  d)  0.020,  e) 0.0020 

1. a)  6dB,  b)  -6dB,  c)  20dB,  d)  -20dB,  e)  40dB. 

2.  a)  240dB,  b)  120dB,  c)  OdB,  d)  -120dB,  e)  -240dB. 
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APPENDIX  B:  INDEX  OF  USEFUL  TABLES  AND  FORMULAS 

Page 

Specific  Acoustic  Impedance  of  Water  &  Air  9 

Pressure  -  P  ^  11 

ref 

Absorption  in  Sea  Water  18 

Speed  of  Sound  in  Sea  Water  22 

Radius  of  Curvature  of  Sound  Ray  26 

Loss  Per  Bounce  in  Surface  Duct  30 

TL  for  Mixed  Layer  31 

Cut  -Off  Frequency  for  the  Mixed  Layer  33 

Convergence  Zone  Range  43 

TL  for  Convergence  Zone  43 

Bottom  Loss  46 

TL  for  Reliable  Acoustical  Path  47 

TL  for  SOFAR  Channel  48 

Band  Levels  from  Pressure  Spectrum  Levels  60 

Nomogram  for  Combining  Levels  63 

Tonals  Combined  with  Background  64 

Typical  Source  Spectrum  Levels  69 

Noise  Spectrum  Level  (Deep  Water)  72 

DI-  for  Line  Array  75 

DI  for  Piston  77 

Typical  ROC  Curve  84 

d  vs  P(FA)  for  P(D)  =  0.5  85 

Energy  Detection  86 

Doppler  Shift  (Passive)  87 

SL  vs  Power  for  Various  DI  95 

Representative  Values  of  TS  98 

Representative  Scattering  Strengths  104 

Doppler  Shift  (Active)  19^ 

Correlation  Detection  19^ 
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